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Foreword 



One of the key features of quantum mechanics is the interference of probability amplitudes. 
The reason for the appearance of interference is mathematically very simple. It is the linear 
structure of the Hilbert space which is used for the description of quantum systems. In terms 
of physics we usually talk about the superposition principle valid for individual and composed 
quantum objects. So, while the source of interference is understandable it leads in fact to 
many counter-intuitive physical phenomena which puzzle physicists for almost hundred years. 

The present thesis studies interference in two seemingly disjoint fields of physics. However, 
both have strong links to quantum information processing and hence are related. In the first 
part we study the intriguing properties of quantum walks. In the second part we analyze a 
sophisticated application of wave packet dynamics in atoms and molecules for factorization 
of integers. 

The main body of the thesis is based on the original contributions listed separately at the 
end of the thesis. The more technical aspects and brief summaries of used methods are left 
for appendices. 
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Part I 

Recurrences in Quantum Walks 
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Introduction 



From random walks to quantum walks 

The term random walk was first introduced by Pearson [Ij in 1905, is a mathematical for- 
malization of a trajectory that consists of successive random steps. Shortly after that a 
paradigmatic application of a random walk - the explanation of Brownian motion [2] and dif- 
fusive processes, was found by Einstein [3] and Smoluchowski [4j. Since then random walks 
have been used in many branches of science [5], ranging from physics, economy, ecology to 
social sciences. Among others, the random walk is one of the cornerstones of theoretical com- 
puter science [61 [7J. Indeed, it can be employed for algorithmic purposes to solve problems 
such as graph connectivity [8], 3-SAT j9] or approximating the permanent of a matrix |10j . 

Quantum walks have been proposed by Aharonov, Davidovich and Zagury [11] as a gener- 
alization of classical random walks to quantum domain. The unitary time evolution governing 
the walk can be either discrete as introduced by Meyer [121 H3] and Watrous [H] leading to 
coined quantum walks or continuous as introduced by Farhi and Gutman [T5| IT6] . It is 
interesting to note that similar ideas can be found already in the works of Feynman [17] 
and Bialynicki-Birula [TS] in the context of discretization of the Dirac equation. Scattering 
quantum walks [T9j [20], |2TJ |22] were proposed by Hillery, Bergou and Feldman as a natural 
generalization of coined quantum walks based on an interferometric analogy. The connection 
between the coined quantum walks and the continuous time quantum walks has been estab- 
lished [2"3"j [21] . Recently, it has been shown that both continuous [23] and discrete time [2B] 
quantum walks can be regarded as a universal computational primitive. By now, quantum 
walks form a well established part of quantum information theory [27]. For a review see e.g. 
the article by Kempe [28] or books by Venegas-Andraca [29J or Konno [30J. 

Continuous-time quantum walks are suitable for the description of coherent transport 
of excitation in networks J2H E2]- Recently, a coherent energy transfer in photosynthetic 
systems was observed [33J. This long-lived coherence which can be described by a gener- 
alized continuous-time quantum walk [34J together with the environmental noise leads to a 
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substantial increase in energy transfer efficiency 

Coined quantum walk is well suited as an algorithmic tool j3U [37] • Several algorithms 
based on coined quantum walks showing speed up over classical algorithms have been pro- 
posed [381 ESI EH HU HI2 H21 HI]. Various properties of coined quantum walks have been 
analyzed, e.g. the effects of the coin and the initial state |1SJ H7J HB], absorbing barriers [15] , 
the hitting times [IHHSOJEI] or the effect of decoherence [131 E2]. Hitting times for continuous 
quantum walks related to the quantum Zeno effect were considered in [53]. Great attention 
has been paid to the asymptotics of quantum walks [5U [571 EH ESI EE]. In particular, local- 
ization was found in 2-D quantum walks [59l 1451 160] and in 1-D for a generalized quantum 
walk [HI1 E2]- Several experimental schemes have been proposed to realize coined quantum 
walks including cavity QED [63] , linear optics [6H [65] , optical lattices [671 [66] , Bose-Einstein 
condensate [68J and quantum rings [69]. Recently, as proof of principle, experiments with 
neutral atoms [70], ions [7T] and photons [72] have been performed. 

In comparison to classical random walks coined quantum walks are considerably more 
flexible. The coin operator can be in principle an arbitrary unitary matrix. Moreover, one 
can choose the initial coin state. All of these influence the dynamics of the quantum walk. 
The diversity of quantum walks asks for a classification. Indeed, in order to exploit the full 
potential of quantum walks for algorithmic purposes one needs to know in which regimes 
they can be operated in. 

The present thesis focuses mainly on one particular quantity which is suitable for the 
classification of both classical as well as quantum walks, namely the probability to return 
to the origin. The recurrence probability is known as the Polya number, after G. Polya 
who as the first discussed this property in the context of classical random walks on infinite 
lattices in 1921 [73J. Polya pointed out the fundamental difference between walks in different 
dimensions. In three or higher dimensions the recurrence probability is less than one and 
depends exclusively on the dimension [74], whereas for walks in one or two dimensions the 
Polya number equals unity. As a consequence, in three and higher dimensions the particle 
has a non-zero probability of escape |75j . Recurrence in classical random walks is closely 
related to first passage times as pointed out in a number of classics papers of statistical 
mechanics [TTJ [76] . A summary of the results on recurrence of classical random walks is left 
for Appendix IA1 

We extend the concept of recurrence and Polya number to quantum walks in Chapter [1] 
based on [I] where a particular measurement scheme was considered. Other possible def- 
initions of the quantum Polya number are briefly discussed following [TTJ . As we show in 
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Appendix [Bj within the framework of our measurement scheme the criterion for recurrence 
of a quantum walk is the same as for the classical random walk - it is determined by the 
asymptotic behaviour of the probability at the origin. To be able to analyze the probability 
at the origin we first solve the time evolution equations. Since the quantum walks in con- 
sideration are translationally invariant we make us of the Fourier transformation and find 
a simple solution in the momentum picture. Probability amplitudes in the position repre- 
sentation are then obtained by performing the inverse Fourier transformation. Hence, they 
have a form of an integral over momenta where the time enters only in the rapidly oscillating 
phase. This allows us to perform the asymptotic analysis of the probability at the origin in 
a straightforward way by means of the method of stationary phase. Basic concepts of this 
method are reviewed in Appendix O We find that the asymptotic scaling of the probability 
at the origin is affected by the additional degrees of freedom offered by quantum mechanics. 
Hence, the recurrence probability of a quantum walk depends in general on the topology of 
the walk, choice of the coin and the initial state. This is in great contrast to classical random 
walks, where the Polya number is characteristic for the given dimension. 

Recurrence of unbiased quantum walks on infinite d- dimensional lattices is analyzed in 
Chapter [2] which is based on [IIIJ. First, we show that for the quantum walk driven by 
Hadamard tensor product coin, the Polya number is independent of the initial conditions, 
thus resembling the property of the classical walks. We provide an estimation of the Polya 
number for this quantum walk in dependence of the dimension of the lattice. Second, we 
examine the Grover walk on a plane, which exhibits localization and thus is recurrent, except 
for a particular initial state for which the walk is transient. We generalize the Grover walk to 
show that one can construct in arbitrary dimensions a quantum walk which is recurrent. This 
is in great contrast with classical random walks which are recurrent only for the dimensions 
d = 1,2. Finally, we analyze the recurrence of the Fourier walk on a plane. This quantum 
walk is recurrent except for a two-dimensional subspace of initial states. We provide an 
estimation of the Polya number in dependence on the initial states. 

In Chapter [3] we extend our analysis of recurrence to biased quantum walks following [IV] . 
As we illustrate in Appendix \A.2\ recurrence of a classical random walk on a line is extremely 
sensitive to the directional symmetry, any deviation from the equal probability to travel in 
each direction results in a change of the character of the walk from recurrent to transient. 
Applying our definition of the Polya number to quantum walks on a line we show that the 
recurrence character of quantum walks is more stable against bias. We determine the range 
of parameters for which biased quantum walks remain recurrent. We find that there exist 
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recurrent genuine biased quantum walks which is a striking difference to classical random 
walks . 

Quantum walks involving more than one particle opens up the possibility of having en- 
tangled initial states or the particles can be indistinguishable - either bosons or fermions. 
In Chapter H] which is based on [Vj we study the motion of two non-interacting quantum 
particles performing a quantum walk on a line. We analyze the meeting problem, i.e. the 
probability that the two particles are detected at a particular position after a certain num- 
ber of steps. The results are compared with the corresponding classical problem which we 
review in Appendix [D] We derive analytical formulas for the meeting probability and find 
its asymptotic behaviour. We show that the decay of the meeting probability is faster than 
in the classical case, but not quadratically as one could expect from the ballistic nature of 
a quantum walk. The effect of non-classical features offered by quantum mechanics on the 
meeting probability is analyzed. We summarize our results and present an outlook in the 
Conclusions. 

Quantum walk on a line - an introductory example 

Before we turn to the presentation of our results we briefly introduce the basic notions of 
quantum walks. For a more comprehensive review we refer to the literature [28] ■ 

Let us begin with the classical random walk on a line. Random walk is a stochastic 
process where the particle moves on an integer lattice in discrete time steps. In each step 
the particle can move from its current location (say m) to the neighboring lattice points (i.e. 
mil) with equal probability. Suppose that the particle is at time t = at the origin of 
the lattice m — 0. After the first step, we can find the particle at site m = 1 or m = -1 
with probability one-half. To calculate the probability that the particle is at position m at a 
latter time t we can use the following recurrence relations 



The solution of the equations ([T]) with the initial condition P(0, 0) = 1 has the form 



Indeed, each random path has the same probability 2~* and the number of paths leading to 
the lattice point m is given by the well-known binomial distribution. It is straightforward to 
calculate various attributes of the random walk, e.g. the mean value and the variance of the 



P(m, t) = -P(m -!,*-!) + -P(m + 1, t - 1) 



m E 1*. 



(1) 




(2) 
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particle's position. We find that the mean value vanishes, in agreement with the unbiasedness 
of the random walk we consider. On the other hand, the variance grows with the square root 
of the number of steps. Indeed, random walk is a diffusion process. 

The quantum walk is a generalization of a classical random walk to a discrete unitary 
evolution of a quantum particle. Hence, there is no randomness in the time evolution itself in 
the quantum case. Nevertheless, the randomness enters through the measurement. Indeed, if 
we want to know the position of the particle we have to measure it and a particular result is 
found with the corresponding probability given by the standard quantum-mechanical formula. 
The particle can be found on any lattice point m G Z. We denote the corresponding position 
eigenstates by \m). These vectors form an orthonormal basis of the position space %p 

Hp = Span {|m)|m e Z} , (m\n) — S mn , \m) (m\ = I. 

m 

As in the classical random walk, the particle moves from its current position to the neigh- 
boring lattice points, but instead of choosing the path randomly it travels all paths simulta- 
neously, i.e. it evolves into a superposition 

\m) — > \m — 1) + \m + 1). 

However, we can easily see that such a time evolution is not unitary. Indeed, two orthogonal 
vectors |0) and |2) evolves into the states 

|0>— >|-1) + |1>, |2>— Hl) + |3>, 

which have non-zero overlap. To make the time-evolution unitary we have to consider a 
particle which has an internal degree of freedom with two orthogonal states \L) and \R). 
This additional degree of freedom is usually referred to as coin and its two orthogonal states 
\L), \R) form a basis of the corresponding coin space %c 

7*c = Span{|L),|i2)}. 

The state of the coin determines the next move of the particle according to 

\m)\L) — ► \m-l)\L), \m}\R) — ► \m+l)\R). 

Such a transformation is performed by the conditional displacement operator S 

S = (\m - l)(m| ® \L)(L\ + \m + l)(m| <g> \R)(R\ ) , 

m 
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which is indeed unitary. However, a time evolution according to S itself would be rather 
trivial. Indeed, if the particle will start the quantum walk in a definite coin state, say \L), 
it will simply move on to the left. Hence, to obtain a non-trivial time evolution we first 
rotate the coin by the coin operator before the conditional displacement S is applied. As 
the coin operator we can in principle choose an arbitrary unitary transformation on the coin 
space He- Here, we consider a particular choice of the Hadamard coin H which performs the 
following rotation 

H \ L ) = -j5 d L > + l*» ' H\R) = ^(\L)~ \R)) ■ 

Finally, we can write the unitary propagator U which performs a single step of the quantum 
walk 

U = S-(I®H). (3) 
Suppose that the particle is initially at the origin with the coin state \L), i.e. 

!#))} = |0}|L). (4) 
After the first step of the quantum walk it evolves into the state 

Ml)> = £/|#))> = ^(|-l)|L) + |l) \R)). (5) 

Note that if we perform the measurement of the particle's position, we find it with equal prob- 
ability on the sites ±1. This is the same result as for the classical random walk. Moreover, 
after the measurement the state of the particle is projected onto the eigenstate corresponding 
to the measurement outcome. Hence, by performing position measurements after each step 
we obtain one classical random path. By making a statistics of such paths we recover a 
classical random walk. To obtain different dynamics we have to let the quantum particle 
evolve unperturbed, i.e. without measurements, for a desired number of steps t, and perform 
the position measurement afterwards. In this way, each path will not obtain probability but 
probability amplitude, which involves a phase. Different paths leading to the same lattice 
point will interfere. Hence, a quantum walk is an interference phenomenon. 

As we have seen from §5§ the probability distribution of the quantum walk after the first 
step does not differ from the probability distribution of the classical random walk. Indeed, if 
the quantum particle is initially localized at the origin no interference can occur. The same 
applies to the second step and the state of the particle is given by 

M2)) = U\m) = l(\- m) + \0)(\L) + \R)) - \2}\R) ) . 
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The probability to find the particle at the position m after two steps P(m, 2) is given by 
P(-2,2) = K-2|(L|^(2))| 2 +|(-2|W(2))| 2 = i 
P(0,2) = |(0|(L|^(2))| 2 + |(0|W(2))| 2 = i 

P(2,2) = |(2|(L|^(2))| 2 + |(2|W(2))| 2 = i 

which is the same as for the classical random walk. Finally, in the third step the interference 
occurs for the first time. The state of the particle after the third step has the form 

iv(3)> = u\m) = ^ (i - m) + 1 - i>(2|l> + ip)) - ii}|p) + is)|p) ) , 

and we see that the probability distribution 

P(-3,3) = |(-3|(L|^(3))| 2 + |(-3|(P|^(3))| 2 = i 
P(-l,3) = |(-1|(L|^(3))| 2 + |(-1|(P|^(3))| 2 = ^ 
P(l,3) = |(1|(L|^(3))| 2 + |(1|(P|^(3))| 2 = ^ 

P(3,3) = |(3|(L|^(3))| 2 + |(3|(P|^(3))| 2 = i 

differs from the classical one. As a consequence of the choice of the initial coin state (jl]) it is 
biased towards the left. 

In general, the state of the particle at a later time t is given by the successive application 
of the propagator U on the initial state |?/>(0)) 

mt)) = u t mo))- (6) 

Let us denote by tpL,(R)(rn,t) the probability amplitude of finding the particle at site m 
with the coin state \L(R)) after t steps of the quantum walk. These amplitudes are the 
coefficients of the decomposition of the state vector \ip(t)) into the basis of the total Hilbert 
space V, = Hp <g> He 

\m) = Yl (Mm,t)\m)\L)+M^,t)\m)\R) ) . (7) 

m 

Using the form of the propagator U ([3]) we find from the time evolution of the state vector 
the equations of motions for the probability amplitudes 

i[) L (m,t) = -^=ij L (m + l,t-l) + -^=ij R (m + l,t-l), 

ip R {m,t) = -^=ip L (m- l,t- 1) - -^=ip R (m- l,t- 1). (8) 
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These equations are reminiscent of the time evolution equations of the classical random 
walk ([I]). However, in (|SJ) we transform probability amplitudes instead of probabilities. The 
probability to find the quantum particle at a particular position m is given by the standard 
quantum-mechanical formula 

P(m,t) = |H(L|^(t)>| 2 + \(m\(Rm))\ 2 = \M^,t)\ 2 + \M™,t)\ 2 . 

In Figure [1] we display the probability distribution of the classical and quantum walk on a 
line obtained from the numerical simulation. Concerning the classical random walk depicted 
by the red points we observe a symmetric gaussian distribution with a rather small width. 
Indeed, the variance of the classical random walk grows with the square root of the number of 
steps, which is a typical signature of diffusion. The probability distribution of the quantum 
walk depicted by the blue points shows striking differences compared to the classical random 
walk. As we have already discussed, due to the choice of the initial coin state the distribution 
is biased to the left. More important observation is that the width of the distribution is 
proportional to the number of steps. Indeed, due to the interference of the probability 
amplitudes (jSj) the growth of the variance is linear in time [54] . Hence, the quantum walk is 
a ballistic process which is the key difference from the diffusive nature of the classical random 
walk. The quadratic speed-up of the variance is at the heart of the fast algorithms based on 
quantum walks [381 ESI SH SH M, S21 Ej • 
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t = 100 




-100 o 100 

m 



Figure 1: The probability distribution of the classical and quantum walk on a line after 100 
steps. For the classical random walk illustrated by the red points we find that the probability 
distribution is peaked at the origin and symmetric. Indeed, the mean value vanishes. The 
width of the distribution is rather small, since the variance of the classical random walk grows 
only with the square root of the number of steps. This is a typical signature of diffusion. 
In contrast, the probability distribution of the quantum walk described by the blue points 
shows striking differences. First, due to the choice of the initial coin state the distribution is 
biased towards left. Second, the width of the distribution is proportional to the number of 
steps. Indeed, the variance of the quantum walk grows linearly with time which is a typical 
signature of a ballistic process. 
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Chapter 1 

Recurrence of Quantum Walks 



Introduction 

Classical random walks are defined as the probabilistic discrete time evolution of the position 
of a point-like particle on a discrete graph. Starting the walker from a well-defined graph 
point (the origin) one can ask whether the particle returns there at least once during the time 
evolution. The probability of this event is called the Polya number [73] . Classical random 
walks are said to be recurrent or transient depending on whether their Polya number equals 
to one, or is less than one, respectively. 

The Polya number of a classical random walk can be defined in the following way 



p = j2^(t), (i.i) 



t=i 

where qo(t) is the probability that the walker returns to the origin for the first time after t 
steps. More practical expression of the Polya number is in terms of the probability po(t) that 
the particle can be found at the origin at any given time instant t. It is straightforward to 
show that 

P = l-— . (1.2) 

EM*) 
t=0 

From (jl.2p we find that the recurrence behaviour of a random walk is determined solely by 
the infinite sum 

oo 

5 = 5>,(t). (1.3) 

t=o 

Indeed, P equals unity if and only if the series S diverges [78]. In such a case the random 
walk is recurrent. On the other hand, if the series S converges, the Polya number P is strictly 
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less than unity and the walk is transient. The well-known result found by Polya [73] is that 
unbiased random walks in one and two dimensions are recurrent while for higher dimensional 
lattices they are transient. For a more detailed review of recurrence of random walks see 
Appendix [A] 

We define the Polya number of a quantum walk in Section 11.11 by considering a specific 
measurement scheme. Other possible measurement schemes are briefly discussed. In accor- 
dance with the classical terminology we describe the quantum walk as recurrent or transient 
depending on the value of the Polya number. We find a condition for the recurrence of a 
quantum walk which is given by the asymptotic behaviour of the probability at the origin. 
A general description of a quantum walk on an infinite d- dimensional lattice is left for Sec- 
tion 11.21 In particular, we find a simple form of the time evolution equation for probability 
amplitudes. In Section 11.31 we employ the translational invariance of the problem which al- 
lows us to solve the equations of motion easily in the momentum representation. We find 
that the amplitudes in the position representation can be written in the form of an integral 
over momenta where the time enters only in the oscillating phase. This form of the solution 
allows a straightforward analysis of the asymptotic behaviour of the amplitudes by means 
of the method of stationary phase. We perform this analysis in Section 11.41 and discuss the 
consequences on the recurrence nature of the quantum walk. In particular, we find that the 
latter is affected by the choice of the coin and the initial coin state. 

1.1 Polya number of a quantum walk 

For quantum walks we can keep the same definition of the Polya number (II. ip being the 
probability of returning to the origin at least once during the time evolution. However, to 
be able to talk about the position of a particle in quantum mechanics one must specify when 
and which type of measurement is performed. According to the definition ( 11.1]) we would 
have to continuously measure whether the particle is at the origin. However, such a radical 
interruption of the system ultimately leads to a loss of coherence which is a vital ingredient 
of a quantum walk. It can be anticipated that within the continuous measurement scheme 
most of the quantum effects become rather weak. The analysis we have performed in [H] 
supports this conclusion. 

In order to preserve the quantum interference we have considered different measurement 
scheme in [I] . The recurrence is understood as a property of an ensemble of particles rather 
than an individual. The measurement scheme is the following: Prepare an ensemble of 
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quantum walk systems in an identical initial state. Take one of such systems, let it evolve for 
one step, perform the measurement at the origin and then discard the system. Take a second, 
identically prepared system, let it evolve for two steps, make a position measurement at the 
origin and then discard it. Continue until a positive outcome is obtained. In the t-th trial 
we do not find the particle at the origin with the probability 1 — po(t). Since the individual 
trials are independent the product 

n 

p n =H(i- P0 (t)) 

t=i 

gives the probability that we have not found any particle at the origin in the first n trials. 
In the complementary event, which occurs with the probability 

n 

P n = l-JI(l-po(t)), (1.4) 

t=i 

we have found at least one particle at the origin. We define the Polya number of a quantum 
walk by extending n to infinity 

+00 

P = l_ JJ(l-po(f)). (1-5) 

t=i 

This definition resembles the expression of the Polya number of a classical random walk in 
terms of the probability at the origin (11. 2ft . However, the inverted sum of po(t) is replaced 
by the product of 1 — Po(t). Nevertheless, we show in Appendix Q3] that definition (ll.5p of the 
Polya number of a quantum walk leads to the same criterion for recurrence in terms of the 
probability at the origin po(t). Indeed, the infinite product in (II. 5p vanishes if and only if the 
series S (11.31) diverges [79] . In such a case the Polya number of a quantum walk is unity and 
we call such quantum walks recurrent. If the series S converges, then the product in (11. 5p 
does not vanish and the Polya number of a quantum walk is less than one. In accordance 
with the classical terminology we call such quantum walks transient. 

The convergence of the series S (11.31) is determined by the asymptotic behaviour of the 
probability at the origin. In the following Sections we find the means which allows us to 
perform this asymptotic analysis. 
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1.2 Description of quantum walks on Z 



Let us first define quantum walks on an infinite d dimensional lattice Z d . The Hilbert space 
of the quantum walk can be written as a tensor product 

U = U P ®U C 

of the position space 

Up = f(Z d ) 

and the coin space Uc- The position space is spanned by the vectors |m) corresponding to 
the particle being at the lattice point m, i.e. 

Up = Span {|m)| m = {mi, . . . , m d } G Z d ] . 

The coin space Uc is determined by the topology of the walk. In particular, its dimension n is 
given by the number of possible displacements in a single step. We denote the displacements 
by vectors 

e, G Z d , i = 1, ... ,n. 

Hence, the particle can move from m to any of the points m + i — 1, . . . , n in a single 
step. We define an orthonormal basis in the coin space by assigning to every displacement 
ej the basis vector je^), i.e. 

U c = Span{|ei)|? = 1, . . . ,n} . 

A single step of the quantum walk is given by 

U = S ■ (I P <g) C) . (1.6) 

Here Ip denotes the unit operator acting on the position space Up. The coin flip operator 
C is applied on the coin state before the displacement S itself. The coin flip C can be in 
general an arbitrary unitary operator acting on the coin space Uc- 

The displacement itself is represented by the conditional step operator S 

S = ^ |m + ej)(m| <g> le^) (e^, 
which moves the particle from the site m to m + e, if the state of the coin is |e^) . 
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Let the initial state of the particle be 

|V(0)> = ^Vi(m,0)|m)<8> 1^}. 

Here ipi(m, 0) is the probability amplitude of finding the particle at time t = at the position 
m in the coin state |e^) . The state of the particle after t steps is given by successive application 
of the time evolution operator given by Eq. ( 11. 6p on the initial state 

|^(t)) = ^Vi(m,t)|m> ® |e,> = E/*|V(0)>. (1.7) 

m,i 

The probability of finding the particle at the position m at time t is given by the summation 
over the coin state, i.e. 

n n 

P ( m ,t) ee ^iHtehKt))! 2 = E l^(m,t)| 2 = t)|| 2 . 

4=1 8=1 

Here we have introduced n-component vectors 

V»(m,t) ee (^ 1 (m,t),V 2 (m,t),...,V„(m,t)) T 

of probability amplitudes. We rewrite the time evolution equation (II. 7p for the state vector 
\if)(t)) into a set of difference equations 

^(m,t) = EC/V'(m-e, J t-l) (1.8) 

for probability amplitudes ?/>(m, 0- Here the matrices C/ have all entries equal to zero except 
for the Z-th row which follows from the coin-flip operator C, i.e. 

(e; \Q\ ej) = 5 u (ei \C\ a,). 

1.3 Time evolution of quantum walks 

The quantum walks we consider are translationally invariant which manifests itself in the fact 
that the matrices C\ on the right-hand side of Eq. (11.81) are independent of m. Hence, the time 
evolution equations (II. 8p simplify considerably with the help of the Fourier transformation 

4>{k,t) = E^(m,t)e im - k , keK d . (1.9) 

m 
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The Fourier transformation defined by Eq. (I1.9P is an isometry between £ 2 (Z d ) and L 2 (K d ) 
where K = (— tt, tt] can be thought of as the phase of a unit circle in R 2 . 

The time evolution in the Fourier picture turns into a single difference equation 

$(k,t) = t7(k)$(k,t-l). (1.10) 

Here we have introduced the propagator in the momentum representation 

f7(k) = D(k) ■ C, D(k) = Diag (e 4ei k , . . . , e 4e "- k ) . (1.11) 

We find that [/ (k) is determined both by the coin C and the topology of the quantum walk 
through the diagonal matrix D(k) containing the displacements e*. 

We solve the difference equation f ll.lOp by formally diagonalising the matrix U(k). Since 
it is a unitary matrix its eigenvalues can be written in the exponential form 

Aj(k) = exp (i ujjik)), 

where the phase is given by the eigenenergy u)j(k). We denote the corresponding eigenvectors 
as Vj(\t). Using this notation the state of the particle in the Fourier picture at time t reads 

tf(M) = E el ^ (%(k),V5(k,0)) Vj (k), (1.12) 

3 

where ( , ) denotes the scalar product in the n dimensional coin space He- Finally, we 
perform the inverse Fourier transformation and find the exact expression for the probability 
amplitudes 

V»(m,t)= / — - d ^(k,t)e— k (1.13) 

in the position representation. 

We are interested in the recurrence nature of quantum walks. As we have discussed in 
Section 11.11 the recurrence of a quantum walk is determined by the asymptotic behaviour of 
the probability at the origin 

Po (t)=p(0,t) = \\^(0,t)f- 

as the number of steps approaches infinity. Hence, we set m = in Eq. (11.131) . Moreover, 
in analogy with the classical problem of Polya we restrict ourselves to quantum walks which 
start at origin. Hence, the initial condition reads 

^(m,0) = 5 m , o ^ V = V>(0,0) (1.14) 
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and its Fourier transformation ip(k, 0) entering Eq. ( I1.12p is identical to the initial state of 
the coin 

^(k,0)=^ 

which is a n-component vector. We note that due to the Kronecker delta in Eq. (11.14p the 
Fourier transformation ip{k, 0) is independent of the momenta k. 

Using the above assumptions we find the exact expression for the probability at the origin 

2 



3=1 



Po(t) 

where Ij(t) are given by the integrals 

/dk 
^ e /,(k), £(k) = ( Vj (k),i;) Vj (k). (1.15) 

K d 

1.4 Asymptotics of the probability at the origin 

Let us discuss how the additional freedom we have at hand for quantum walks influences 
the asymptotics of the probability at the origin po(t). We suppose that the functions u)j(k) 
and fj(k) entering Ij(t) are smooth. According to the method of stationary phase [80] which 
we briefly review in Appendix O the major contribution to the integral Ij(t) comes from 
the stationary points k° of the eigenenergies u)j(k), i.e. from the points where the gradient 
vanishes 

Vw 7 -(k) = 0. 
J k=k° 

The asymptotic behaviour of Ij (t) is then determined by the stationary point with the greatest 
degeneracy given by the dimension of the kernel of the Hessian matrix 

dk m dk n 

evaluated at the stationary point, i.e. by the flatness of Uj(k). The function fj(k) entering the 
integral I jit) determines only the pre-factor. We now discuss how the existence, configuration 
and number of stationary points affect the asymptotic behaviour of Ij(t). As a rule of thumb, 
the decay of the probability at the origin po{t) can slow down with the increase in the number 
of stationary points. Let us briefly discuss the results. 
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1.4.1 No stationary points 

If ujj(k) has no stationary points then Ij(t) decays faster than any inverse polynomial in t. 
Consequently, the decay of the probability at the origin is also exponential 

Po(t) ~ e-* 

with some positive rate 7. Indeed, quantum walks for which the probability at the origin 
decays so fast are transient. Such a situation occurs e.g. for extremely biased quantum walks 
which we analyze in Chapter [3j 

1.4.2 Finite number of stationary points 

Suppose that tOj(k) has a finite number of non-degenerate stationary points, i.e. the deter- 
minant of the Hessian matrix H is non-zero for all stationary points. If the function /j-(k) 
does not vanish at the stationary points then the contribution from all stationary points to 
the integral Ij(t) is of the order t~ d l 2 . Consequently, the probability at the origin behave like 

Po (t) ~ r d 

as t approaches infinity. Clearly, the sum S defined in (11.31) is convergent for d > 1. Hence, 
the quantum walks for which the eigenenergies have only non-degenerate stationary points 
are recurrent only for the dimension d — 1, i.e. on a line. This is e.g. the case of the 
Hadamard walk with tensor product coin studied in Chapter 12.11 

1.4.3 Continuum of stationary points 

If 0Jj(k) has a continuum of stationary points then the dimension of the continuum determines 
the decay of the integral Ij(t). The case of 2-D integrals with curves of stationary points are 
treated in [SO]- It is shown that the contribution from the continuum of stationary points to 
the integral Ij(t) is of the order t~ l l 2 . This is greater than the contribution arising from a 
discrete stationary point which is of the order t~ l . Hence, the continuum of stationary points 
has effectively slowed-down the decay of the integral Ij(t). Consequently, the leading order 
term of the probability at the origin is 

Po(t) -r 1 , 

and we find that such a quantum walk is recurrent. We come across this situation in the case 
of the Fourier walk on a plane in Chapter 12.31 Similar results can be expected for higher 
dimensional quantum walks where u>j(k) have a continuum of stationary points. 



24 



A special case for a continuum of stationary points is when Wj(k) does not depend on 
n variables, say ki, . . . k n , but has a finite number of stationary points with respect to the 
remaining d — n variables k n+ i, . . . ,kd- Indeed, such an Uj(k n +i, ■ ■ ■ ,kd) has obviously a zero 
derivative with respect to ki, % = 1, . . .n. Suppose that the function fj(k) factorizes 



fj( k ) = 9j(h, ■ ■ ■ , k n ) ■ hj(k n+1 , ...,h 



d 



In such a case Ij(t) is given by the product of time- independent and time-dependent integrals 
over n and d — n variables 



' /k -e l ^ k ^-^%{k n+ll ... 1 k d ) 



(2-ny-n 

.K d ~ n 

It is easy to find that if the time-independent integral does not vanish Ij(t) behaves asymp- 
totically like t~( d ~ n ^ 2 . Hence, the asymptotic behaviour of the probability at the origin 
is 

The quantum walks of this kind would be recurrent if the eigenenergy Uj would depend only 
on a single component of the momenta k. In the extreme case when cjj(k) does not depend 
on k at all we can extract the time dependence out of the integral Ij(t). If the remaining time 
independent integral does not vanish then po(t) converges to a non-zero value and say that 
such a quantum walk exhibits localization. Note that since po(t) has a non- vanishing limit the 
quantum walk is recurrent. Indeed, localization implies recurrence. We find localization in 
Chapter 12.21 for the Grover walk on a plane. Moreover, extending the 2-D Grover walk to Z d 
we find quantum walks where some of the eigenenergies are either constant or depend only 
on a single momentum component. As discussed above, such quantum walks are recurrent. 

1.4.4 Effect of the initial state 

So far we have assumed that the function /j(k) is non- vanishing for k values corresponding 
to the stationary points. However, the initial state ip can be orthogonal to the eigenvector 
Vj{\s) for k = k° corresponding to the stationary point. In such a case the function /j(k) 
vanishes for k = k° and the stationary point k° does not contribute to the integral Ij{t). 
Consequently, the decay of po(t) can speed up. Hence, for quantum walks we might change 
the recurrence behaviour and the actual value of the Polya number by altering the initial 
state ip. Indeed, we find this non-trivial effect of the initial state for the Grover walk and the 
Fourier walk on a plane in Chapters 12.21 and 12.31 
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Chapter 2 



Recurrence of Unbiased Quantum 
Walks on Infinite Lattices 

Introduction 

In the present Chapter we determine the recurrence behaviour and the Polya number of 
several unbiased quantum walks. We concentrate on the effect of the coin operators and the 
initial states. For this purpose we fix the topology of the walks. We consider quantum walks 
where the displacements have all entries equal to ±1 

ei = (1, • • • , 1) T , • • • , e 2 d = (-1, . . . , -1) T . 

In such a case the coin space has the dimension n = 2 d where d is the dimension of the 
lattice. Moreover, the diagonal matrix D(k) entering the propagator in the Fourier picture 
(11. lip can be written as a tensor product 

D(k) = D(h) <g> . . . ® D(k d ) (2.1) 

of 2 x 2 diagonal matrices 

D(kj) = Diag (> ' 7 ' <V 7 ' ) . 

This fact allows us to extend some of the results for the quantum walks on a line or on a 
plane to quantum walks on a d- dimensional lattice. 

First, in Section I2TT1 we treat Hadamard walk on Z d with an independent coin for each 
spatial dimension. We find that for this quantum walk the probability at the origin is 
independent of the initial coin state. Hence, a unique Polya number can be assigned to 
this quantum walk for each dimension d. In contrast with the classical random walks the 
Hadamard walk is recurrent only for d = 1. In Section 12.21 we analyze the recurrence of 
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the Grover walk on a plane. This quantum walk exhibits localization [HD] and therefore is 
recurrent. However, for a particular initial state localization disappears and the Grover walk 
becomes transient. We find an approximation of the Polya number for this particular initial 
state. We then employ the Grover walk on a plane to construct for arbitrary dimension d a 
quantum walk which is recurrent. This is in great contrast with the classical random walks, 
which are recurrent only for the dimensions d — 1,2. Finally, in Section I2~3l we analyze the 
Fourier walk on a plane. This quantum walk is recurrent except for a two-parameter family 
of initial states for which it is transient. For the latter case we find an approximation of the 
Polya number depending on the parameters of the initial state. We summarize our results 
in Section 



2.1 Hadamard walk on 



Let us start with the analysis of the recurrence behaviour of the Hadamard walk on a line 
which is driven by the coin 

We find that the propagator in the Fourier picture 

U H {k) = D{k) ■ H = _l-4k j (2.2) 

has eigenvalues e % where the phases u>i(k) are given by 

ru\ ■ f sink \ fu \ . {sink 

oji(k) = arcsm I — — , cu 2 (/c) = — 7r — arcsm 



Thus the derivatives of Ui with respect to k reads 

dui(k) diij2{k) cos A; 

dk dk ^2 - sin 2 k 



(2.3) 



and we find that the phases Ui(k) have common non-degenerate stationary points k° = ±7r/2. 
It follows that the probability at the origin behaves asymptotically like This asymptotic 
scaling is independent of the initial state. Indeed, no non-zero initial state ip exists which 
is orthogonal to both eigenvectors at the common stationary points k° = ±7r/2. Hence, the 
Hadamard walk on a line is recurrent, i.e. the Polya number equals one, independent of the 
initial coin state. 
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We turn to the Hadamard walk on a d- dimensional lattice. The coin flip operator has the 
form of the tensor product of d 2 x 2 Hadamard matrices 



H d = H g> ...®H. (2.4) 

Hence, we have an independent coin for each spatial dimension. It follows that also the 
propagator in the Fourier picture has the form of the tensor product 

U Hd (k) = U H (h) <g> . . . ® U H (k d ) (2.5) 

of d time evolution operators given by Eq. ( 12. 2 p with different momenta k{. Hence, the 
eigenenergies of the propagator (12. 5p have the form of the sum 

d 

Wi (k) = J>;^)- ( 2 - 6 ) 

i=i 

Therefore we find that the asymptotic behaviour of this quantum walk follows directly from 
the asymptotics of the Hadamard walk on a line. Indeed, the derivative of the phase w 3 -(k) 
with respect to ki reads 

(2.7) 



duj(k) duj t (ki) 



dki dki ' 

and so w?(k) has a stationary point k° = (k®, k®, ■ ■ ■ , k d ) if and only if for all / = 1, . . . , d the 
point /c z ° is the stationary point of ojj^ki, oci). As we have found from Eq. (12.31) the stationary 
points of Uj t are kf = ±7r/2. Hence, all phases ujj(k) have 2 d common stationary points 
k° = (±7r/2, . . . , ±7r/2). It follows that the asymptotic behaviour of the probability po(t) is 
given by 

Po(t) ~ t~ d . (2.8) 

As follows from the results for the Hadamard walk on a line the asymptotic behaviour given 
by Eq. (12. 8p is independent of the initial coin state. Compared to classical walks this is a 
quadratically faster decay of the probability at the origin which is due to the quadratically 
faster spreading of the probability distribution of the quantum walk. 

We illustrate the results for Hadamard walk on a plane driven by the coin 

/ 1 1 1 1 \ 
1-1 1-1 
1 1-1-1 

\i -l -l 1/ 

in Figure 12.11 Here we show the probability distribution in dependence on the initial state 
and the probability at the origin po(t). The first two plots indicates that the initial state of 



H 2 = i 
2 



(2.9) 
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the coin influences mainly the edges of the probability distribution. However, the probability 
Po(t) is unaffected and is exactly the same for all initial states. The lower plot confirms the 
asymptotic behaviour of the probability at the origin po(t) ~ t~ 2 . 

Since the probability at the origin po(t) decays like t~ d we find that the Hadamard walk 
on Z d is recurrent only for dimension d — 1 and is transient for all higher dimensions d > 2. 
Moreover, the whole sequence of probabilities Po(t) is independent of the initial state. Hence, 
the Polya number for this class of quantum walks depends only on the dimension of the walk 
d, thus resembling the property of the classical walks. On the other hand, this quantum 
walk is transient for the dimension d = 2 and higher. This is a direct consequence of the 
faster decay of the probability at the origin which, in this case, cannot be compensated for 
by interference. 

Let us estimate the value of the Polya number for the dimension d > 2. As depicted 
in the lowest plot of Figure 12.11 the probability at the origin approaches quite rapidly its 
asymptotic form 

Mt) M -jjj. (2.10) 

Hence, already the first few terms of the product in Eq. (11. 4p are sufficient to estimate the 
value of the Polya number. Taking into account the first three terms of po(t) which are found 
to be 

Pd(2) = i p (4)=po(6) = ^, (2.11) 
we obtain the following approximation of the Polya number 




(2.12) 



We compare the estimation in Eq. (I2.12p with the numerical results obtained from the 
simulation of the Hadamard walk with 1000 steps in the Table 12.11 and find that they are in 
excellent agreement. 



2.2 Grover walk on a plane 



We turn to the Grover walk on a plane which is driven by the coin 



G 



1 

2 



/ -1 
1 
1 
1 



1 \ 
1 
1 

1 / 



(2.13) 
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t 



Figure 2.1: Probability distribution of the Hadamard walk on a plane after 50 steps and the 
probability at the origin po(t) for different choices of the initial state. In the upper plot we 
choose the initial state |(1, i, i, — 1) T which leads to a symmetric probability distribution, 
whereas in the middle plot we choose the initial state (1,0, 0,0) T resulting in a dominant 
peak of the probability distribution in the lower-left corner of the (m, n) plane. However, 
the initial state influences the probability distribution only near the edges. The probability 
Po(t) is unaffected and is the same for all initial coin states. The lower plot confirms the 
asymptotic behaviour of the probability at the origin po(t) ~ t~ 2 independent of the initial 
state. 
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Dimension 


Simulation 


Estimation (]2.l2|) 


Error (%) 


2 


0.29325 


0.27325 


6.8 


3 


0.12947 


0.12841 


0.82 


4 


0.06302 


0.06296 


0.01 


5 


0.031313 


0.031309 


0.01 



Table 2.1: Comparison of the Polya number for the Hadamard walk on Z d obtained from 
the numerical simulation and the estimation of Eq. (12.1 2p . 



It was identified numerically [45] and later proven analytically [60] that the Grover walk 
exhibits a localization effect, i.e. the probability po(t) does not vanish but converges to a 
non-zero value except for a particular initial state 

^ G = Vg(0, 0, 0) = l - (1, -1, -1, 1) T . (2.14) 

In order to explain the localization we analyze the eigenvalues of the propagator in the 
Fourier picture for the Grover walk 

U G (k l ,k 2 ) = {D(k 1 )®D(k 2 ))G. (2.15) 

We find that they are given by 

A lj2 = ±l, \ 3A (k 1 ,k 2 ) = e ±i ^ k ^ (2.16) 

where the phase u(ki, k 2 ) reads 

cos(a;(A;i, k 2 )) = — cos k\ cos k 2 . (2-17) 

The eigenvalues Ai,2 are constant. As a consequence the probability at the origin is non- 
vanishing as discussed in detail in Chapter ll.4.3[ unless the initial state is orthogonal to 
the eigenvectors corresponding to Xi 2 at every point (ki,k 2 ). By explicitly calculating the 
eigenvectors of the matrix Uc{ki, k 2 ) it is straightforward to see that such a vector is unique 
and equals that in Eq. (I2.14p . in agreement with the result derived in |60j . 

It is easy to show that for the particular initial state given by Eq. (I2.14p the probability 
Po(t) decays like t~ 2 . Indeed, as the initial state of Eq. (I2.14p is orthogonal to the eigenvectors 
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corresponding to Ai^ the asymptotic behaviour is determined by the remaining eigenvalues 
A 3i4 (fci, k 2 ), or more precisely by the stationary points of u(ki, k 2 ). From Eq. ( 12.1 7p we find 
that it has only non- degenerate stationary points k®, k\ = ±7r/2. For the initial state of 
Eq. f)2.14p the probability that the Grover walk returns to the origin decays like t~ 2 . We 
conclude that the Grover walk on a 2-D lattice is recurrent and its Polya number equals one 
for all initial states except the one given in Eq. (I2.14p for which the walk is transient. We 
illustrate these results in Figure 12.21 and Figure 12.31 

In Figure 12.21 we show the probability distribution generated by the Grover walk and the 
probability at the origin for a symmetric initial state 

^s=\{l,i,i,-l) T . (2.18) 

This particular choice of the initial state results to a probability distribution with a dominant 
central spike, as depicted in the upper plot. The lower plot indicates that the probability at 
the origin has a non- vanishing limit. 

In contrast for the initial state ipc given by (12.141) the central spike in the probability 
distribution vanishes, as we illustrate in the upper plot of Figure [2731 The lower plot indicates 
that the probability at the origin decays like t~ 2 . 

Let us estimate the Polya number of the Grover walk for the initial state of Eq. (12.141) . 
The numerical simulations indicate that the probability at the origin po(t) for the initial state 
ipa is t ne same as the probability at the origin of the 2-D Hadamard walk. Hence, their Polya 
numbers coincide. With the help of the relation (I2.12p we can estimate the Polya number of 
the Grover walk with the initial state of ipc by 

P G ^ G ) = P H2 « 0.27325. (2.19) 

The above derived results allow us to construct a quantum walk which is recurrent for 
an arbitrary dimension d, except for a subspace of initial states. Let us first consider the 
case when the dimension of the walk is even and equals 2d. We choose the coin as a tensor 
product 

G 2d = ® d G (2.20) 

of d Grover coins given by Eq. (12.1 3p . As follows from Eqs. (11.111) and (12.11) the time 
evolution operator in the Fourier picture is also a tensor product 

U G2d (k) = U G (h, k 2 ) ® . . . ® U G (hd-i, k 2d ) (2.21) 
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Figure 2.2: Probability distribution of the Grover walk after 50 steps and the probability at 
the origin for a symmetric initial state (I2.18p . This particular choice of the initial state leads 
to a symmetric probability distribution with a dominant central spike, as depicted in the 
upper plot. The lower plot indicates that the probability at the origin has a non- vanishing 
limit as t approaches infinity. The results are qualitatively the same for all initial coin states 
except for ipc given in (I2.14p , as we illustrate in Figure 12.31 
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Figure 2.3: Probability distribution of the Grover walk after 50 steps and the probability 
at the origin for a particular initial state ipG given by Eq. (I2.14p . In contrast to Figure [2721 
we find that the central spike vanishes and most of the probability is situated at the edges. 
Moreover, the probability at the origin vanishes as t approaches infinity, as we illustrate in 
the lower figure. Here we plot the probability po(t) multiplied by t 2 to unravel the asymptotic 
behavior of the probability at the origin. The plot confirms the analytic result of the scaling 

P o(t) ~ r 2 . 



34 



of the matrices JJq defined by Eq. ( 12.1 5p with different Fourier variables hi. Hence, the 
eigenvalues of Ua 2d (k) are given by the product of the eigenvalues of Uq- Since two eigenvalues 
of Ug are constant as we have found in Eq. (12 . 1 6[) one half of the eigenvalues of £/G 2d (k) are 
also independent of k. As we have discussed in Chapter 11.4.31 the probability po(t) converges 
to a non-zero value and therefore the quantum walk exhibits localization. 

In the case of odd dimension 2c? + 1 we augment the coin given by Eq. (I2.2(jp by the 
Hadamard coin for the extra spatial dimension 

GWi = G 2d <g> H. (2.22) 

Performing a similar analysis as in the case of even dimensions we find that for the quantum 
walk driven by the coin G2d+i the probability that the walk returns to the origin decays like 
t -1 due to the Hadamard walk in the extra spatial dimension. Hence, this quantum walk is 
recurrent. 

We note that due the fact that the 2-D Grover walk is transient for the initial state ipa 
the same statement holds for the above constructed quantum walks, supposed that the initial 
state contains ipc in its tensor product decomposition. Such vectors form a subspace with 
dimension equal to 4 d_1 for even dimensional walks and 2 x 4 d_1 for odd dimensional walks. 

2.3 Fourier walk on a plane 

We turn to the 2-D Fourier walk driven by the coin 



(2.23) 



As we will see, the Fourier walk does not exhibit localization. However, the decay of the 
probability po(t) is slowed down to t^ 1 so the Fourier walk is recurrent, except for a subspace 
of states. 

We start our analysis of the Fourier walk with the propagator 

U F (k 1 ,k 2 ) = (D(k 1 )®D(k 2 ))F, 

which determines the time evolution in the Fourier picture. It seems to be hard to determine 
the eigenvalues of U F (ki, k 2 ) analytically. However, we only need to determine the stationary 
points of their phases Uj(ki, k 2 ). For this purpose we consider the eigenvalue equation 

k 2 , u) = det (u F (k u k 2 ) - e i u 'ij = 0. 
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This equation gives us the eigenenergies 0Ji{k\, k 2 ) as the solutions of the implicit function 



k 2 , oS) — 1 + cos(2A; 2 ) — 2 cos(2o;) + 2 sin 2a; + 4 cos k 2 sin a; (sin k\ — cos fci) = 0. 

Using the implicit differentiation we find the derivatives of the phase uj 
du cos k 2 si n w ( cos &i + srn k\ ) 



du 
dko 



cos(2w) + sin(2o;) + cos k 2 cos a; (sin &i — cos fc 2 ) 

2 sin fc 2 sinw (cos k\ — sin fci) — sin(2&; 2 ) 
2 (cos(2a;) + sin(2a;) + cos k 2 cosw (sin k\ — cos fc 2 )) 



(2.24) 



with respect to k\ and k 2 . Though we cannot eliminate u on the RHS of Eq. (I2.24p . we can 
identify the stationary points k° = (k®, k 2 ) 



du (k) 



dh 



0. 



1,2 



k=k° 



of o;(fci, fe) with the help of the implicit function $(&i, &2, w). We find the following: 
(i) Wi i2 (fci,/c 2 ) have stationary lines 



7i = (fci,0) and 7 2 = (&i,7r) 



(m) all four phases Ui{k\, k 2 ) have stationary points for 

fc° = -, - — and A: 2 ° = ±- 
1 4 4 2 2 

It follows from the discussion of Chapter 11.4.31 that the two phases £^1,2(^1,^2) with sta- 
tionary lines 7i j2 are responsible for the slow down of the decay of the probability po(t) to t^ 1 
for the Fourier walk, unless the initial coin state is orthogonal to the corresponding eigenvec- 
tors v 1,2(^1, k 2 ) at the stationary lines. For such an initial state the probability po(t) behaves 
like t~ 2 as the asymptotics of the integral given by Eq. fl 1 . 1 5 [) is determined only by the 
stationary points (ii). 

Let us determine the states ipp which lead to the fast decay t~ 2 of the probability that 
the Fourier walk returns to the origin. The states ipF have to be constant vectors fulfilling 
the conditions 

(1*1,2 (k), ip F ) = V k e 71,2, 
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which implies that ipp must be a linear combination of ify^k G 7^2) forming a two-dimensional 
subspace in % c . For k 2 = 0,n we can find the eigenvectors of the matrix Up(ki, k 2 ) explicitly 

if 



Ul(Ki,UJ 


= ^2(^1, 


IT) 




ui(fci,7r) 


= ^2(^1, 


0) 


= i(-e- lfc M,e-* 


u 3 (fci,0) 


= "3(^1, 


7T) 


= -^(1,0,1, Of 


^4(^1,0) 


= ^4(^1, 


7T) 


= -1(0, 1,0, -if. 


The explicit form of if) f reads 












(a, 6, a, — &f , 



(2.25) 



where a, 6 6 C. We point out that the particular initial state 



/ 1, 1-A 1 / 1 — i 1-A 

which was identified in [45J as the state which leads to a symmetric probability distribution 
with no peak in the neighborhood of the origin belongs to the family described by Eq. f)2.25p . 

We illustrate the results in Figure 12.41 and Figure 12.51 In Figure 12.41 we plot the prob- 
ability distribution and the probability po(t) for the Fourier walk with the initial state 
if} = (1, 0, 0, 0) T . This vector is not a member of the family if>F(a, b) defined by Eq. f)2.25p . 
We find that a central peak is present, as depicted in Figure |2~41 However, in contrast to the 
Grover walk, the peak vanishes as shown by plotting the probability po(t) multiplied by t in 
Figure I2~41 indicating a decay like t~ l , in agreement with the analytical result. In contrast, 
for Figure 12.51 we have chosen the initial state given by Eq. ( I2.26P which is a member of the 
family -0j?(a,6). The upper plot shows highly symmetric probability distribution. However, 
the central peak is not present and as the lower plot indicates the probability po(t) decays 
like r 2 . 

We conclude that the Fourier walk is recurrent except for the two-dimensional subspace 
of initial states defined by Eq. (12.251) for which the walk is transient. 

We turn to the estimation of the Polya numbers of the 2-D Fourier walk for the two- 
dimensional subspace of initial states given by Eq. (12.251) . We make use of the normalization 
condition and the fact that the global phase of a state is irrelevant. Hence, we can choose a 
to be non-negative real and b is then given by the relation 



a?e H 
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Figure 2.4: Probability distribution after 50 steps and the time evolution of the probability 
Po(t) for the Fourier walk with the initial state if) = (1,0,0,0) T . The upper plot of the 
probability distribution reveals a presence of the central peak. Indeed, ip is not a member 
of the family ipp{a,b). However, in contrast to the Grover walk the peak vanishes. In the 
lower plot we illustrate this by showing the probability Po(t) multiplied by t to unravel the 
asymptotic behaviour po(t) ~ t" 1 . 
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Figure 2.5: Probability distribution after 50 steps and the time evolution of the probability 
Po(t) for the Fourier walk with the initial state given by Eq. f !2.26p . Since ip is a member of 
the family i/jf(ci, b) the central peak in the probability distribution is not present, as depicted 
on the upper plot. The lower plot indicates that the probability po(t) decays like t~ 2 . 
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Therefore, we parameterize the family of states defined by Eq. (I2.25P by two real parameters 
- a ranging from to ^ and the mutual phase G [0,2-71"). The exact expression for 
Po{a, 0, t) can be written in the form 

Ki(t) — K 2 {t)aJ^ ~ a 2 (cos0 — sin0) 
p (a,(f>,t) = - , 

where K\ 2 has to be determined numerically. Nevertheless, the numerical simulation of 
Po(a,(f),t) at two values of (a, 0) enables us to find the numerical values of Kx t2 (t) and we 
can evaluate po{a, 0, t) at any point (a, 0). The probability po(a, 0, t) shows the maximum at 
a — ~, (j) — and the minimum for the same value of a and the phase = ^f. Consequently, 
these points also represent the maximum and the minimum of the Polya numbers. 

In Figure 12.61 we present the approximation of the Polya number Eq. f ll.4p in its depen- 
dence on a and and a cut through the plot at the value a = 1/2 containing both the global 
minimum and the global maximum. Here we have evaluated the first 100 terms of po{a, 0, t) 
exactly. We see that the values of the Polya number vary from the minimum p™ n pa 0.314 
to the maximal value of pp ax pa 0.671. We note that for the initial states that do not belong 
to the subspace defined by Eq. (I2.25P the Polya number equals one. 

2.4 Conclusions 

Our results, summarized in Table I2.2[ demonstrate that there is a remarkable freedom for 
the value of the Polya number for quantum walks, depending both on the initial state and 
the coin operator, in contrast to the classical random walk where the dimension of the lattice 
uniquely defines the recurrence probability. Hence, the quantum Polya number is able to 
indicate physically different regimes in which a quantum walk can be operated in. 
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Figure 2.6: Approximation of the Polya numbers for the 2-D Fourier walk and the initial 
states from the family of states defined by Eq. (12.251) in their dependence on the parameters 
of the initial state a and <fi. Here we have evaluated the first 100 terms of Po(a, <fr, t) exactly. 
The Polya numbers cover the whole interval between the minimal value of P™ n « 0.314 
and the maximal value of pp ax pa 0.671. The extreme values are attained for a = 1/2 and 
cf) mm = 7n/4, respectively <fi max = 3n/4. On the lower plot we show the cut at the value 
a = 1/2 containing both the maximum and the minimum. 
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Quantum walk 


Section 


Probability at the origin 


Polya number 


d-D Hadamard 




j.~d ^ 


1 for d = 1 
< 1 for d > 2 
independent of if) 


2-D Grover 


|2.2| 


const, for ifj ^ ipc 
t~ 2 for if) = if)Q 


1 

< 1 


2-D Fourier 


12.31 


t^ 1 for if) j£ if>p 
t~ 2 for if) Eif>F 


1 

< 1 

dependent on if) 



Table 2.2: Summary of the main results. We list the types of studied quantum walks, the 
asymptotic behaviour of the probability at the origin and the Polya number in the respective 
cases in its dependence on the initial state if). 
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Chapter 3 

Recurrence of Biased Quantum Walks 
on a Line 

Introduction 

Recurrence of classical random walks is a consequence of the walk's symmetry. As we briefly 
review in the Appendix I A. 2 \ they are recurrent if and only if the mean value of the position of 
the particle vanishes. This is due to the fact that the spreading of the probability distribution 
of the position is diffusive while the mean value of the position propagates with a constant 
velocity. In contrast, for quantum walks both the spreading of the probability distribution 
and the propagation of the mean value are ballistic. In the present Chapter we show that 
this allows for maintaining recurrence even when the symmetry is broken. 

The Chapter is organized as follows: In Section I3TT1 we describe the biased quantum walk 
on a line, find the propagator in the momentum representation and solve the time evolution 
equation. The recurrence of the quantum walk is determined by the asymptotics of the 
probability at the origin. We perform this analysis in Section 13.21 and find the conditions 
under which the biased quantum walk on a line is recurrent. In Section 13.31 we analyze 
the recurrence of biased quantum walks from a different perspective. We find that the 
recurrence is related to the velocities of the peaks of the probability distribution of the 
quantum walk. The explicit form of the velocities leads us to the same condition derived in 
Section 13.21 Finally, in Section 13.41 we derive the formula for the mean value of the position 
of the particle in dependence of the parameters of the walk and the initial state. We find that 
there exist genuine biased quantum walks which are recurrent. We summarize our results in 
the conclusions of Section 13.51 
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3.1 Description of the walk 



Let us consider biased quantum walks on a line where the particle has two possibilities - 
jump to the right or to the left. Without loss of generality we restrict ourselves to biased 
quantum walks where the jump to the right is of the length r and the jump to the left has a 
unit size. We depict the biased quantum walk schematically in Figure 13.11 



\L) \R) 




-2-10 1 r r + 1 



Figure 3.1: Schematics of the biased quantum walk on a line. If the coin is in the state \R) 
the particle moves to the right to a point at distance r. With the coin state \L) the particle 
makes a unit length step to the left. Before the step itself the coin state is rotated according 
to the coin operator C(p). 



The Hilbert space of the particle has the form of the tensor product 

U = U P ® U c 

of the position space 

Hp = f(Z d ) = Span{|m)| m G Z} , 
and the two dimensional coin space 

U C = C 2 = Span{|i?),|L)}. 

The propagator of the quantum walk in the position representation is 

U = S (J P <g> C) , 

where the displacement operator S has the form 

+oo +oo 

S = ^2 \m + r}(m\®\R)(R\+ ^ \m-l)(m\®\L)(L\. 

m=— oo m=— oo 

The coin flip C can be in general an arbitrary unitary operator acting on the coin space 
He- However, as has been discussed in [IS] the probability distribution is not affected by 
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the complex phases of the coin operator. Hence, it is sufficient to consider the one-parameter 
family of coins 

From now on we restrict ourselves to this family of coins. The value of p = 1/2 corresponds 
to the well known case of the Hadamard walk. 

In the momentum representation the propagator has the form 

U(k) = Diag (e* e-) . C(p) = ( J^_, t ) . 

Since it is a unitary operator its eigenvalues are e l UJ1 ' 2 where the phases are given by 

Ui(k) = — - — k + arcsin ^/psm kjj , 

r — 1 / /r + l\\ 

U2{k) = — - — k — it — arcsin I i/psin I — - — k . (3.1) 



2 

We denote the corresponding eigenvectors by V\^{k). We give their explicit form in Sec- 



tion 13.51 The solution of the time evolution equation in the Fourier picture has the standard 
form 

2 

?(M) = E ei u ' m (*i(*)>^(M)) vj(k), 

where ip(k, 0) is the Fourier transformation of the initial state. We restrict ourselves to the 
situation where the particle is initially localized at the origin as dictated by the nature of 
the problem we wish to study. Hence, the Fourier transformation of such an initial condition 
is equal to the initial state of the coin which we denote by ip. Since ip can be an arbitrary 
normalized complex two-component vector we parameterize it by two parameters a G [0, 1] 
and if G [0, 2it) in the form 



^(v»> (3 - 2) 

The solution in the position representation is obtained by performing the inverse Fourier 
transformation 

1>(rn,t) = f § £(M) e~ mk = £ f (vj(k),ip)vj(k). (3.3) 
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3.2 Asymptotics of the probability at the origin 



To determine the recurrence nature of the biased quantum walk we have to analyze the 
asymptotic behaviour of the probability at the origin. Exploiting (13. 3p the amplitude at the 
origin reads 

mt) = Y,j f^ <w '«* (Vjikl^vjik), (3.4) 

which allows us to find the asymptotics of the probability at the origin with the help of the 
method of stationary phase [80]. The important contributions to the integrals in (13. 4p arise 
from the stationary points of the phases ( I3.ip . We find that the derivatives of the phases 
LO\^{k) are 

u/( k ) = r ~ l I VEfr + 1) cos (fcffi 
1 2 v /4 + 2p[cos(Jfe(r + l))-l]' 

u/(k) = ^ + l)cos(fc^) (35) 

2 ^4 + 2p [cos(k(r + 1)) - 1] ' 

Using the method of stationary phase we find that the amplitude will decay slowly - like t~^, 
if at least one of the phases has a vanishing derivative inside the integration domain. Solving 
the equations u[ 2 {k) = we find that the possible stationary points are 



2 / (l-p)(r-l) 2 \ , N 

k = ± arccos ± J ^ . 3.6 

r + 1 \ V P r I 

The stationary points are real valued provided the argument of the arcus-cosine in ( 13. 6 j) is 
less or equal to unity 

Apr 

This inequality leads us to the condition for the biased quantum walk on a line to be recurrent 

TlY- (3-7) 
r + 1 J 

We illustrate this result in Figure I3T21 for a particular choice of the walk parameter r = 3. 

Our simple result proves that there is an intimate nontrivial link between the length of 
the step of the walk and the bias of the coin. The parameter of the coin p has to be at least 
equal to a factor determined by the size of the step to the right r for the walk to be recurrent. 
We note that the recurrence nature of the biased quantum walk on a line is determined only 
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h , , , J 

-7T -f f n 

k 

Figure 3.2: The existence of stationary points of the phases coi^k) in dependence on the 
parameter p and a fixed step length r. We plot the implicit functions u[ 2 (k) = for r = 3. 
The plot indicates that for p < pn(3) = \ the phases ui\j,{k) do not have any stationary 
points. Consequently, the probability amplitude at the origin decays fast and such biased 
quantum walk on a line is transient. For p > Pr(3) the stationary points exists and the 
quantum walk is recurrent. 

by the parameters of walk itself, i.e. the coin and the step, not by the initial conditions. The 
parameters of the initial state a and ip have no effect on the rate of decay of the probability 
at the origin. 

3.3 Velocities of the peaks 

We can determine the recurrence nature of the biased quantum walk on a line from a dif- 
ferent point of view. This approach is based on the following observation. The well known 
shape of the probability distribution generated by the quantum walk consists of two counter- 
propagating peaks. In between the two dominant peaks the probability is roughly indepen- 
dent of m and decays like t~ l . On the other hand, outside the decay is exponential as we 
depart from the peaks. As it has been found in [51] the positions of the peaks varies linearly 
with the number of steps. Hence, the peaks propagate with constant velocities, say vl and 
Vr. For the biased quantum walk to be recurrent the origin of the walk has to remain in 
between the two peaks for all times. In other words, the biased quantum walk on a line is 
recurrent if and only if the velocity of the left peak is negative and the velocity of the right 
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peak is positive. 

The velocities of the left and right peak are easily determined. We rewrite the formula 
(13. 3p for the probability amplitude ip(m,t) into the form 

iP(m,t) = J2j Y /^ {k) - ak)t (v,(k)^)v 3 (k), 

where we have introduced a = y. Due to the fact that we concentrate on the amplitudes at 
the positions m ~ t we have to consider modified phases 

ujj(k) = ojj{k) — ak. 

The peak occurs at such a position mo where both the first and the second derivatives of 
uij(k) vanishes. The velocity of the peak is thus ao = n f-. Hence, solving the equations 

= r-l Vp( r + l)cos(^) 

1 2 ^A + 2p [cos(fc(r + 1)) - 1] 
~, m r-l y/p(r + 1) cos (ffi) 

W 9 (k) = — Q = 0, 

2 2 ^4 + 2p [cos(fc(r + 1)) - 1] 

y/2[2- p + pcos{k{r + l)))2 

for a determines the velocities of the left and right peak Vl,r- The third equation is inde- 
pendent of a and we easily find the solution 

Ann . 2n(2n + 1) 

k = — — , k = — , n E Z. 

r + 1 r + 1 

Inserting this ko into the first two equations we find the velocities of the left and right peak 

r-l r + r-l r + l_ /ooN 

v l = — g Y~ v^' ^ = + ~~ 2~ v 7 ^- ( 3 - 8 ) 

We illustrate this result in Figure 13.31 where we show the probability distribution generated 
by the quantum walk for the particular choice of the parameters r = 3, p = The initial 
state was chosen according to a = ^ and tp = ir. Since the velocity of the left peak vl is 
negative this biased quantum walk is recurrent. 

The peak velocities have two contributions. One is identical and independent of p, the 
second is a product of r and p and differs in sign for the two velocities. The obtained results 
indicate that biasing the walk by having the size of the step to the right equal to r results 
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rn 

Figure 3.3: Velocities of the left and right peak of the probability distribution generated by 
the biased quantum walk on a line and the recurrence. We have chosen the parameters r = 3, 
a = p = -t= and <p = tt. The left peak propagates with the velocity vl ~ —0.68, the velocity 
of the right peak is vr 2.68. In between the two peaks the probability distribution behaves 
like t^ 1 while outside the decay is exponential. Since the velocity vl is negative the origin 
of the walk remains in between the left and right peak. Consequently, this quantum walk is 
recurrent. 

in dragging the whole probability distribution towards the direction of the larger step. This 
is manifested by the term £=- which appears in both velocities with the same sign. On 
the other hand the parameter of the coin p does not bias the walk. As we can see from the 
second terms entering the velocities it rather influences the rate at which the walk spreads. 

As we have discussed above the biased quantum walk on a line is recurrent if and only if 
Vl is negative and vr is positive. The form of the velocities (I3.8P implies that this condition 
is satisfied if and only if the criterion (13. 7p is fulfilled. 

3.4 Mean value of the position 

As we discuss in the Appendix IA.2I the classical random walks are recurrent if and only if 
the mean value of the position vanishes. We show that this is not true for biased quantum 
walks, i.e. there exist biased quantum walks on a line which are recurrent but cannot produce 
probability distribution with zero mean value. This is another unique feature of quantum 
walks compared to the classical ones. 
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Let us derive the formula for the mean value of the position of the particle for the biased 
quantum walk. With the help of the weak limit theorem [HS] we express the mean value after 
t steps in the form 



• i */ — 7T 



t 

up to the corrections of the order 0(t~ l ). Here Vj(k) are eigenvectors of the unitary propa- 
gator U(k), u)j(k) are the derivatives of the eigenenergies and ifi is the initial state expressed 
in (13. 2p . The derivatives of the phases are given in (I3.5p . We express the eigenvectors in the 
form 

(A;) = n 1 (A ; )(VW,-Vp + e i(wiW - rfc) ) T , 

v 2 (k) = n 2 (k) (yr^,-Vp + e^ 2W - r - fc) ) T . 

The normalizations are given by 

ni{u) = 2 — 2^/pcos (u — arcsin [y/psinu]) , 
n 2 (u) = 2 + 2y/~p cos {u 4~ arcsin [y/p sin u\) , 

where we have introduced u = Ml±i) ^ Q shorten the notation. The mean value is thus given 
by the following integral 



x\ f f(a,<p,p,r,u)du + 0(r u 

t/ J 2(r + 1)-7T [l + y^coswx] [l — v / r psinu 2 \ 





where 

u\ = u + arcsin (y/p sin u), u 2 = u + arccos( v / psin-u), 
and the numerator reads 

f(a,ip,p,r,u) = (1 — p) [r — 1 + p(a + r(a — 1)) (1 + cos(2*u)) + 



+ y/a(l — a)y/p(l - p)(r + 1) (cosy? + cos(y? + 2m)) 
Performing the integrations we arrive at the following formula for the position mean value 

(f) ~ (I" V / W)(«(r + !)-!) + ^ V / W + 

y^a(l - o)(l - VI - p)(l - p)(r + 1) cosy , 
H ). (d.yj 

VPl 1 -p) 
50 



We see that for quantum walks the mean value is affected by both the fundamental walk 
parameters through r and p and the initial state parameters a and tp. The mean value is 
typically non- vanishing even for unbiased quantum walks ( with r = 1 ). However, one 
easily finds [45J that the initial state with the parameters a = 1/2 and if = tt/2 results in 
a symmetric probability distribution with zero mean independent of the coin parameter p. 
Indeed, the quantum walks with r = 1, i.e. with equal steps to the right and left, does not 
intrinsically distinguish left from right. On the other hand the quantum walks with r > 1 
treat the left and right direction in a different way. Nevertheless, one can always find for 
a given r a coin parameter p such that for all p > po the quantum walk can produce a 
probability distribution with zero mean value. This is impossible for quantum walks with 
p < Po and we will call such quantum walks genuine biased. 

Let us determine the minimal value of p for a given r for which mean value vanishes. We 
first find the parameters of the initial state a and <p which minimizes the mean value. Clearly 
the term on the second line in (13. 9p reaches the minimal value for ip = it. Differentiating 
the resulting expression with respect to a and setting the derivative equal to zero gives us 
the condition 

2 | (2a-l)y/^T r P) _ 
Py/ a{\ — a) 

on the minimal mean value with respect to a. This relation is satisfied for Og = |(1 — yfp)- 
The resulting formula for the mean value reads 

x\ r-1 t (l-y/l-p-p){l +_r) 



1 1 aom 2 2^J (1 — p)p 

This expression vanishes for 

/r 2 -l\ 2 

Po(r) = [—^ • (3.H) 

Since (13. lOj) is a decreasing function of p the mean value is always positive for p < p 
independent of the choice of the initial state. For p > po one can achieve zero mean value for 
different combination of the parameters a and (p. 

The formula (13.111) is reminiscent of the condition (j3.7p for the biased quantum walk on a 
line to be recurrent. However, r is in (13. 1 ip replaced by r 2 . Therefore we find the inequality 
Pr < Po- Hence, the quantum walks with the coin parameter Pr < p < po are recurrent 
but cannot produce a probability distribution with zero mean value. We conclude that there 
are genuine biased quantum walks which are recurrent in contrast to situations found for 
classical walks. 
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3.5 Conclusions 



We have analyzed one dimensional biased quantum walks. Classically, the bias leading to a 
non-zero mean value of the particle's position can be introduced in two ways — unequal step 
lengths or unfair coin. In contrast, for quantum walks on a line the initial state can introduce 
bias for any coin. On the other hand, for symmetric initial state modifying only the unitary 
coin operator while keeping the equal step lengths will not introduce bias. Finally, the bias 
due to unequal step lengths may be compensated for by the choice of the coin operator for 
some initial conditions. For this reason we have introduced the concept of the genuinely 
biased quantum walk for which there does not exists any initial state leading to vanishing 
mean value of the position. 

We have determined the conditions under which one dimensional biased quantum walks 
are recurrent. This together with the condition of being genuinely biased give rise to three 
different regions in the parameter space which we depict as a " phase diagram" in Figure 13.41 




r 

Figure 3.4: "Phase diagram" of biased quantum walks on a line. The horizontal axis repre- 
sents the length of the step to the right r and the vertical axis shows the coin parameter p. 
The dotted line corresponds to the recurrence criterion (13.71) . while the squares represent the 
condition (13. lip on the zero mean value of the particle's position. The quantum walks in the 
white area are transient and genuine biased. In between the two curves (light gray area) we 
find quantum walks which are recurrent but still genuine biased. The quantum walks in the 
dark gray area are recurrent and for a particular choice of the initial state they can produce 
probability distribution with vanishing mean value. 
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Chapter 4 

Meeting Problem in the Quantum 
Walk 

Introduction 

In this Chapter we study the evolution of two particles performing a quantum walk. The 
evolution of each of the two particles is subjected to the same rules. One of the interesting 
questions, when two particles are involved, is to clarify how the probability of the particles 
to meet changes with time (or number of steps taken in walk). Because the behavior of a 
single particle performing a quantum walk differs from its classical counterpart it has to be 
expected that the same applies to the situation when two particles are involved. Interference, 
responsible for the unusual behavior of the single particle should play also a considerable role 
when two particles are involved. The possibility to change the input states (in particular 
the possibility to choose entangled initial coin states) adds another interesting point to the 
analysis. In the following, we study the evolution of the meeting probability for two particles. 
We point out the differences to the classical case and discuss the influence of the input state 
on this probability. 

Before we turn to the meeting problem we generalize in Sections 14. II and 14.21 the quantum 
walk to two distinguishable and indistinguishable particles. The meeting problem for two 
distinguishable particles with initially factorized coin states is analyzed in Section 14.31 We 
derive the asymptotic behavior of the meeting probability and compare it with the results for 
the classical random walk which are summarized in Appendix |D] The effect of entanglement 
on the meeting probability is considered in Section 14.41 Finally, in Section 14.51 we analyze 
the meeting probability for two indistinguishable bosons and fermions. We summarize our 
results in the conclusions of Section 14.61 
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4.1 Quantum walk with two distinguishable particles 



The Hilbert space of the two particles is given by a tensor product of the single particle 
spaces, i.e. 

H = (Hp® H c )i ® (Hp <g> Hc^- 
Each particle has its own coin which determines his movement on the line. Since we assume 
that there is no interaction between the two particles they evolve independently and the 
time evolution of the whole system is given by a tensor product of the single particle time 
evolution operators. We describe the state of the system by vectors 



i[>(m, n, t) 



( ipLL(m,n,t) \ 
TpR L {m,n,t) 
TpLii{m,n,t) 



where e.g. the component ipRL(m,n,t) is the amplitude of the state where the first particle 
is on m with the internal state \R) and the second particle is on n with the internal state 
\L). The state of the two particles at time t is then given by 

1^(0) = 5^ ^(^M^^MHk j'h- (4.1) 

m,n i,j= R,L 

The conditional probability that the first particle is on a site m at time t, provided that the 
second particle is at the same time at site n, is defined by 

P(m,n,t)= |<m,i|<n,j|^(t)>| 2 = ^ |^(m, n, t)\ 2 . (4.2) 

i,j=L,R i,j=L,R 

Note that if we would consider a single quantum particle but on a two dimensional lat- 
tice, with two independent Hadamard coins for each spatial dimension, ( 14. 2p will give the 
probability distribution generated by such a two dimensional walk. This shows the relation 
between a one dimensional walk with two particles and a two dimensional walk. The reduced 
probabilities for the first and the second particle are given by 

P 1 (m, t) = J2 P K n, t), P 2 (n, t) = ^ P(m, n, t). 

n m 

The dynamics of the two particles is determined by the single particle motion. Since we 
can always decompose the initial state of the two particles into a linear combination of a 
tensor product of a single particle states and because the time evolution is also given by a 
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tensor product of two unitary operators, the shape of the state will remain unchanged. Thus 
we can fully describe the time evolution of the two quantum particles with the help of the 
single particle wave-functions. A similar relation holds for the probability distribution (14. 2 j) . 
Moreover, in the particular case when the two particles are initially in a factorized state 

|V»(0)> = ( XJVu(™,0)|m,i>i J ® ( ^^(n,0)|n,j) 2 J , (4.3) 

\m,i / \ nj / 

which translates into ipij(m,n,0) = i/;ii(m,0)ip2j(n,0). Hence, the probability distribution 
remains a product of a single particle probability distributions 

P{m,n,t) = (|^iL(m,t)| 2 + |V li? (m,t)| 2 )(|V 2 L(n,t)| 2 + |V2i?(n,t)| 2 ) 

= P 1 (m,t)P 2 (n,t). (4.4) 

However, when the initial state of the two particles is entangled 

l#))> = Ei (E^( m '°)l m '*M ® E*°)KjM L (4.5) 

the probability distribution cannot be expressed in terms of single particle distributions but 
probability amplitudes 

2 



P(m,n,t) = ^ 

i,j=L,R 



X>5(tM)V$(M) 



(4.6) 



Notice that the correlations are present also in the classical random walk with two parti- 
cles, if we consider initial conditions of the following form 

P(m,n,0) = ^Pr(m,0)P 2 >,0). (4.7) 

a 

The difference between ( 14. 6 p and (14. 7p is that in the quantum case we have probability am- 
plitudes not probabilities. The effect of the quantum mechanical dynamics is the interference 
in 1Q]1 . 

Let us define the meeting problem. We ask for the probability that the two particles will 
be detected at the position m after t steps. This probability is given by the norm of the 
vector ifj(m, m, t) 

M D {m,t)= I^J 

(m,m,t)| 2 = P(m,m,t). (4.8) 

i,j=L,R 
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As we have seen above for a particular case when the two particles are initially in a factorized 
state of the form H4.3|) this can be further simplified to the multiple of the probabilities that 
the individual particles will reach the site. However, this is not possible in the situation when 
the particles are initially entangled (14. 5p . The entanglement introduced in the initial state 
of the particles leads to the correlations between the particles position and thus the meeting 
probability is no longer a product of the single particle probabilities. 



4.2 Quantum walk with two indistinguishable particles 

We analyze the situation when the two particles are indistinguishable. Because we work with 
indistinguishable particles we use the Fock space and creation operators, we use symbols at ~ 
for bosons and b}„ for fermions, e.g. at creates one bosonic particle at position m with 
the internal state The dynamics of the quantum walk with indistinguishable particles is 
defined on a one-particle level, i.e. a single step is given by the following transformation of 
the creation operators 

(m,L) [K \ (m— 1,L) (m+l,R) I ' (m,R) [K \ (m—l,L) (m+l,R) I ' 



^2 V (m-l,L) 1 "(m+l.fl) J ' "(m.-R) ' 

for bosonic particles, similarly for fermions. The difference is that the bosonic operators 
fulfill the commutation relations 







while the fermionic operators satisfy the anticommutation relations 

b(m,i),b( n j) }=0, |o(m,i),o[ nJ) | = 5 mn 5ij. (4.10) 

We will describe the state of the system by the same vectors of amplitudes ip(m,n,t) as for 
the distinguishable particles. The state of the two bosons and fermions analogous to ( 14.1 j) 
for two distinguishable particles is given by 

I^b(*)> = ^ m ^ n ^) h \m,i) h \n,j)\ VaC )^ 
m,n i,j=L,R 

\$F{t)) = Yl Yl ^j( m > n > t )Hm4)Hn,j)\ VaC )' ( 4 - 1X ) 
m,n i,j=L,R 

where \vac) is the vacuum state. Note that in (14.111) both summation indexes m and n 
run over all possible sites, even though e.g. the vectors at -.a) Jvac) and at ^at .Avac) 
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correspond to the same physical state. Using the commutation ( 14. 9 p and anticommutation 
(I4.10p relations we can restrict the sums in ( 14. lip over an ordered pair (m, n) with m > n. 
The resulting wave-function will be symmetric or antisymmetric. 
The conditional probability distribution is given by 

P BtF {m,n,t)= ^ \( 1 (m,i) 1 (n,j)\^B,F{ t ))\ 2 = IV'ij (^, n, t) ± ^ («, m, t) | 2 , 

i,j=L,R i,j=L,R 

for m 7^ n, and for m = n 

I 1 2 i 1 2 i 1 2 

P B (m,m,t) = |(2 (OT)L) |V' j b(*))| + \(2{ m , R )\ip B {t))\ + | (l( m ,L)l( m ,R)\ipB(t))\ 

= 2 \^ LL (m,m,t)\ 2 + 2 \ip RR (m,m,t)\ 2 + \ip LR (m,m,t) + tp RL (m,m,t)\ 2 
= M B {m,t), 

I 1 2 2 

P F (m,m,t) = \(l( m ,L)l(m,R)WF(t))\ = \fpLR(m,m,t) - ip RL (m,m,t)\ 

= M F (m,t). (4.12) 

The diagonal terms of the probability distribution (I4.12p define the meeting probability we 
wish to analyze. 

Let us specify the meeting probability for the case when the probability amplitudes can 
be written in a factorized form ^ (m, n, t) = ipnim, t)ip2j{n, t), which for the distinguishable 
particles corresponds to the situation when they are initially not correlated. In this case the 
meeting probabilities are given by 

M B (m, t) = 2 |Vii(m, t)^ 2L (m, t)\ 2 + 2 \ip 1R (m, t)ip 2R (m, t)\ 2 

+ \ip 1L (m,t)i) 2 R(m,t) + ip 1R (m,t)ip 2 L(m,t)\ 2 , (4.13) 

for bosons and 

M F (m,t) = \ipiL{m,t)ip2R{m,t) - ?p 1R (m,t)tJj2L(m,t)\ 2 , (4.14) 

for fermions. We see that they differ from the formulas for the distinguishable particles, 
except for a particular case when the two bosons start in the same state, i. e. ipi{ m iO) = 
ip2{m-i 0) = ijj(m, 0) for all integers m. For this initial state we obtain 

M B (m } t) = |^(m,t)| 4 + |^(m,t)| 4 + 2|^ L (m,t)^(m,t)| 2 
= (|^(m,t)| 2 + |^(m,t)| 2 ) 2 

which is the same as for the case of distinguishable particles starting at the same point with 
the same internal state. 
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4.3 Meeting problem for distinguishable particles 



Let us compare the meeting problem in the classical and quantum case. We study the two 
following probabilities: the total meeting probability after t step have been performed 

M(t) = ^M(m,t), (4.15) 

m 

and the overall meeting probability during some period of steps T defined as 

T 

M(T) = 1- JJ(1-M(*)). (4.16) 

t=i 

The total meeting probability M(t) is the probability that the two particles meet at time t 
anywhere on the lattice, the overall meeting probability A4(T) is the probability that they 
meet at least once anywhere on the lattice during the first T steps. 

We first concentrate on the influence of the initial state on the meeting probability for the 
distinguishable particles. We consider three situations, the particles start localized with some 
initial distance 2d (for odd initial distance they can never meet, without loss of generality 
we assume that the first starts at the position zero and the second at the position 2d), with 
the coin states: 

(i) right for the first particle and left for the second 

fe(0,2d,0) = 1, 

(ii) symmetric initial conditions l/y/2(\L) +i\R)) for both 

#),2d,0) = ^ \ , 

V- 1 / 

(Hi) left for the first particle and right for the second 

^ Lfl (0,2d,0) = 1. 

In the first case the probability distributions of the particles are biased to the right for the 
first particle, respectively to the left for the second, and thus the particles are moving towards 
each other. In the second case the particles mean positions remain unchanged, as for this 
initial condition the single particle probability distribution is symmetric and unbiased. In 
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the last case the particles are moving away from each other as their probability distributions 
are biased to the left for the first one and to the right for the second. 

Let us specify the meeting probabilities (14.151) . Since the two particles are initially in a 
factorized state it follows from (I4.4p and (14. 8 p that the meeting probability is fully determined 
by the single particle probability distribution. Let 

l^ (L) W) = ^(^WjK^+tfWjK-R)) (4.17) 

m 

|^)(t)) = ^(^f(m,t)|m,L)+^ ) (m,t)|m, J R)) (4.18) 

m 

be the state of a single quantum particle after t steps, under the assumption that the initial 
condition was 

\^ L \0)) = \0,L), \^ R \0)) = \0,R). 

Let us denote by p( L ' R \m, t) the corresponding single particle probability distributions. The 
meeting probabilities for the three situations (i)-(iii) are then given by 

M RL {t,d) = ^P (iJ) (m,t)P (i) (m-2d,t) 
fit 

^ P< L ) (m, t) + pW (m, t) P< L ) (m - 2d, t) + P (L) (m - 2d, t) 
M s {t, d) = 2 2 

m 

M LR {t,d) = ^P^(m,t)P (Ji) (m-2d,t). (4.19) 

m 

Figure 14.11 shows the time evolution of the meeting probability for the three studied 
situations and compares it with the classical case. The initial distance is set to and 10 
lattice points. The plot clearly shows the difference between the quantum and the classical 
case. 

In contrast to the classical walk, in the quantum case the meeting probability is oscillatory. 
The oscillations arise from the single particle probability distribution. After some rapid 
oscillations in the beginning we get a periodic function with the characteristic period of 
about six steps, independent of the initial state. In the quantum case the maximum of the 
meeting probability is reached sooner than in the classical case - the number of steps needed 
to hit the maximum is linear in the initial distance d. This can be understood from the shape 
of the particles probability distribution. The maximum of the meeting probability is obtained 
when the peaks of the probability distribution of the first and second particle overlap. If the 
initial distance between the two particles is 2d then the peaks will overlap approximately 
after \/2d steps. The value of the maximum depends on the choice of the initial state. 
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Initial distance d = 





Figure 4.1: Time evolution of the meeting probability for the three types of initial states 
and the classical random walk with two particles. The initial distance is set to (upper 
plot) and 10 lattice points (lower plot). The upper plot shows a faster decay of the meeting 
probability when the two particles are initially at the same lattice point. Indeed, quantum 
walk spreads quadratically faster compared to the classical random walk. Since both particles 
start the walk from the origin the results for the initial states \LR) and \RL) are identical. 
In the lower plot, where the particles are initially separated, we observe an increase in the 
meeting probability for quantum walk. On the other hand, on a long-time scale the meeting 
probability decays faster in the quantum case. 
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We turn to the meeting probabilities on a long-time scale. We present the review of the 
results derived in Appendix [D] For the classical random walk we find in Appendix ID.ll that 
the meeting probability can by estimated by 

M d (t, d)^^= exp(-^) ~ -±=(1 - (4.20) 



for large number of steps t. We see that the asymptotic behaviour of the meeting probability 
is determined by t~^. Concerning the quantum walk, in Appendix ID. 21 we approximate the 
single particle probability distribution according to [54J and replace the sum in (14.191) by 
integral. We find that within this approximation the meeting probability can be expressed in 
terms of the elliptic integrals. Finally, using the asymptotic expansion of the elliptic integrals 
we find the behaviour of the meeting probability for large number of steps 

In 

M D (t,d) (4.21) 

Hence, the meeting probability decays faster in the quantum case compared to the classical 
case f)4.20p . However, the decay is not quadratically faster, as one could expect from the fact 
that the single particle probability distribution spreads quadratically faster in the quantum 
walk. The peaks in the probability distribution of the quantum walk slow down the decay. 

Note that the estimation (14.211) holds for d > 0, i.e. the initial distance has to be non- 
zero. As we mention in Appendix ID.2[ the continuous approximation of the single particle 
probability distribution is not quadratically integrable, and therefore we cannot use this 
approach for the estimation of the meeting probability when the two particles are initially at 
the same lattice point. There does not seem to be an easy analytic approach to the problem. 
However, from the numerical results, the estimation 

M D (t) ~ ^ (4.22) 

fits the data the best. 



We illustrate these results in Figure 14.21 We plot the meeting probability multiplied by 
the number of steps to unravel the different scaling in the classical and quantum case. In the 
upper plot both particles start from the origin, whereas in the lower plot the initial distance 
is 10 lattice points. The numerical results are consistent with the analytical estimation of 
(14.211) and support the approximation (14. 22p . 

We focus on the overall meeting probability defined by (14.161) . In Figure 14.31 we plot the 
overall probability that the two particles will meet during the first T = 100 steps. 
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Figure 4.2: Long-time behaviour of the meeting probability in the classical and quantum 
walk. In the upper plot both particles start the walk from the origin. In the lower plot the 
particles are initially separated by 10 lattice points. To highlight the asymptotic scaling of 
the meeting probability we plot the latter one multiplied by the number of steps. We can 
clearly see the difference between the classical and quantum walk. In the quantum case the 
re-scaled meeting probability shows a logarithmic increase. On the other hand, the growth is 
much faster (with a square root of t) for the classical case. The numerical results are in good 
agreement with the analytical estimation of Appendix |D] which are summarized in (I4.2ip . 
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Figure 4.3: The overall meeting probability for two distinguishable quantum and classical 
particle during first 100 steps as a function of the initial distance. The same plot on the 
logarithmic scale. Only the values for even points are plotted since for odd initial distance 
the particles never meet. 
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On the first plot we present the difference between the three studied quantum situations, 
whereas the second plot, where the meeting probability is on the log scale, uncovers the 
difference between the quantum and the classical random walk. In the log scale plot we see 
that the overall meeting probability decays slower in the quantum case then in the classical 
case, up to to the initial distance of \/2T. This can be understood by the shape and the 
time evolution of a single particle probability distribution. After t steps the maximums of 
the probability distribution are around the point s ± where s is the initial starting point 
of the quantum particle. For t = 100 steps the peaks are around the points s ± 70. When 
the two particles are initially more then 140 points away, the peaks do not overlap, and the 
meeting probability is given by just the tails of the single particle distributions, which have 
almost classical behavior. 

Finally, we note that the overall meeting probability Ai (T, d) defined in (14.161) converges 
to one as T approaches infinity for both classical and quantum walk independent of the initial 
distance. Indeed, to estimate hA (T, d) we rewrite it in the form 



M (T, d) = 1 -exp 



lnm(l-M(M)) 



v t=l 



(4.23) 



and estimate the exponent with the first order Taylor expansion 

/ T \ T T 

In [11(1- M(t, d)) = J2 ln (! " M ^ d )) ~ - E M & ( 424 ) 



,t=i / t=i t=i 



The scaling of the meeting probability M(t,d) both in the classical case (I4.20p and in the 
quantum case (14.21 j) is slow enough such that the sum in (I4.24p diverges to — oo as T grows. 
Consequently, the exponential in (I4.23P vanishes as T grows. Hence, the overall meeting 
probability converges to unity for both classical and quantum walk, i.e. the particles will 
meet with certainty during their time evolution. 



4.4 Effect of the entanglement 

We will consider the case when the two distinguishable particles are initially entangled. Ac- 
cording to (14.61) the meeting probability is no longer given by a product of a single particle 
probability distributions. However, it can be described using single particle probability am- 
plitudes. We consider the initial state of the following form 

|V>(0)) = |0,2d)®|x>, 
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where \x) is one of the Bell states 

l^> = -±=(\LR)±\RL)) : 

|0±) = -L(\LL)±\RR)). (4.25) 
The corresponding probability distributions resulting from such initial states have the form 
P^(m,n,t) = ~ \4 L \™,t)i;f\n-2d,t)±4 R \m,t)^ L \n-2d,t)\\ 

i,j=L,R 

P^(m,n,t) = - \^ L \m,t)^\n-2d,t)±4 R \m,t)tpf\n-2d,t) 2 ,(4.26) 

i,j=L,R 

where ijj L ^ R \m,t) are the probability amplitudes from (I4.18P which describe the state of a 
single particle after t steps starting the quantum walk from the origin with the initial coin 
state L(R). The meeting probabilities are given by the sum of the diagonal terms in (I4.26P 

■ 2 

o V V m L '(m,. fMb (R) (m, - 2d. t) + ib (R) (m,. tU\ L '(m. - 2d. t) 



M^±(t,d) = \J2 \4 L \rnMf\rn-2d,t)±4 R \m,t)^ L \m-2d,t) 

m i,j=L,R 

M^±(t,d) = \4 L \m,t)^ L \m-2d,t)±4 R \m,t)^ R \m-2d,t) 



m i,j=L,R 

The reduced density operators for both coins are maximally mixed for all four Bell states 
f)4.25p . From this fact follows that the reduced density operators of the particles are 

Pi® = \(\¥ L \t))(^ L \t)\ + \^ R \t)){^ R \t)\) 

= l(\4 L \t))(4 L \t)\ + \4 R \t))(4 B) (t)\ 

where the states it)) are analogous to \ib L ( R \t)) expressed in (14.181) but with shifted 

starting point by 2c?, i.e. 

\4 L ' R \t)) = (^L' R \rn - 2d,t)\m,L)+^ R \m - 2d,t)\m,R) 

m 

The reduced probabilities are therefore 

Pi(m,t) = ^(P (L) (m,t) + P {R) (m,t)) 

P 2 (m,t) = P 1 {m-2d,t), (4.27) 

which are symmetric and unbiased. Notice that the product of the reduced probabilities 
(14.271) gives the probability distribution of a symmetric case studied in the previous section. 



65 



Therefore to catch the interference effect in the meeting problem we compare the quantum 
walks with entangled coin states (I4.25P with the symmetric case Ms- Figure 14.41 shows 
the meeting probabilities and the difference M x — Ms, the initial distance between the two 
particles was chosen to be 10 points. 

We see that the effect of the entanglement could be both positive or negative. Notice 
that 

M^-(t,d)-M s (t,d) = -{M (t>+ {t,d)-M s {t,d)) 
M^-(t,d)-M s {t,d) = -{M^+{t,d)-M s {t,d)), 

so the effect of is opposite to and \<f>~) is opposite to The main difference is 

around the point t w \/2d, i.e., the point where for the factorized states the maximum of the 
meeting probability is reached. The peak value is nearly doubled for M^-, but significantly 
reduced for M^+. On the long time scale, however, the meeting probability M^- decays faster 
than in the other situations. According to the numerical results presented in Figure I4.5[ the 
meeting probabilities for and \<f> ) maintain the asymptotic behavior lnt/t, but for 
it goes like 

M^-(t,d) ~ -. 

The initial entanglement between the particles influences the height of the peaks giving the 
maximum meeting probability and affects also the meeting probability on the long time scale. 

Let us briefly comment on the overall meeting probability. As we have discussed in the 
previous section the overall meeting probability converges to one only if the decay of the 
meeting probability is not faster than |. As we have seen the entanglement could speed-up 
the decay of the meeting probability but it is never faster than |. Therefore we conclude 
that for the initially entangled particles the overall meeting probability converges to one. 

4.5 Meeting problem for indistinguishable particles 

We turn to the meeting problem for two indistinguishable particles. As an example, we 
consider the initial state of the form |l(o,R)l(2d,L)), i-e. one particle starts at the site zero 
with the right coin state and one starts at 2d with the left state. This corresponds to the 
case Mrl for the distinguishable particles. The meeting probabilities are according to (I4.13p . 
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Initial distance d = 10 




Figure 4.4: Comparison of the meeting probability for the initially entangled coins and the 
symmetric case. The initial distance between the two particles is set to 10 points. As the 
initial coin states we choose the Bell states f)4.25p . We observe that the effect of the entangled 
coin state on the meeting probability can be both positive or negative. In the lower plot we 
show the difference in the meeting probability with respect to the symmetric case. We find 
that the effect of is opposite to and is opposite to 
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Initial distance d = 10 
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t 

Figure 4.5: Asymptotic behaviour of the meeting probability for the initially entangled coins. 
In order to unravel the asymptotic scaling of the meeting probability we multiply M x (t) by 
the number of steps t. We see that for the Bell states \ip + ) (black dots) and l^) (open 
circles/diamonds) the rescaled meeting probability M x (t) ■ t shows a logarithmic increase 
with t, while for (stars) the value of M x (t) ■ t levels. These results indicate that the 
asymptotic decay of the meeting probability is faster for the singlet state compared to 
the other Bell states or factorized initial conditions. 

04.141) given by 

M B (t,d) = ]T(2|#W)| 2 |4 L V-2^ 

m 

+| (m, t)^f (m - 2d, t) + if,™ (m, f ) V#° ( m ~ ^ *) ? ) . 
M F (t,d) = ^(|^ ) (m,t)^ ) (m-2d,t)-#(m,t)^ L) (m-2d,t)| 2 ^ (4.28) 

m 

In Figure fl~6l we plot the meeting probabilities and the difference Mb,f — Mrl- 
From the figure we infer that the peak value is in this case only slightly changed. Signifi- 
cant differences appear on the long time scale. The meeting probability is greater for bosons 
and smaller for fermions compared to the case of distinguishable particles. This behavior can 
be understood by examining the asymptotic properties of the expressions (14.281) . Numerical 
evidence presented in Figure 14.71 indicates that the meeting probability for bosons has the 
asymptotic behavior of the form ln(t)/£. However, for fermions the decay of the meeting 
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Initial distance d = 10 





Figure 4.6: Comparison of the meeting probability for bosons, fermions and distinguishable 
particles. The initial distance between the two particles is set to 10 points. We find that the 
maximum value of the meeting probability is almost unaffected. However, for longer times 
we observe an increase in the meeting probability for bosons and decrease for fermions. In 
the lower plot we show the difference in the meeting probability for bosons and fermions with 
respect to distinguishable particles. We find that the increase of the meeting probability for 
bosons is the same as the decrease for fermions. 
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probability is faster having the form 

M F (t,d) ~ j. 

The fermion exclusion principle simply works against an enhancement of the meeting prob- 
ability. 



Initial distance d = 10 




t 

Figure 4.7: Asymptotic behaviour of the meeting probability for bosons, fermions and distin- 
guishable particles. In order to unravel the asymptotic scaling of the meeting probability we 
multiply M(t) by the number of steps t. We see that for bosons (stars) and distinguishable 
particles (black dots) the rescaled meeting probability M(t) ■ t shows a logarithmic increase 
with t, while for fermions (open circles) the value of M(t) ■ t levels. These results indicate 
that the meeting probability decays faster for fermions. 

For the overall meeting probability we can use the same arguments as in the previous 
section and conclude that it will converge to one for both bosons and fermions. 

4.6 Conclusions 

We have defined and analyzed the meeting problem in the quantum walk on an infinite line 
with two quantum particles. For distinguishable particles we have derived analytical formulas 
for the meeting probability. The asymptotic behavior following from these results shows that 
the meeting probability decays faster but not quadratically faster than in the classical random 
walk. This results in the slower convergency of the overall meeting probability, however it 
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still converges to one. This is due to the fact that the meeting probability does not decay 
faster than |. Such a situation might occur in higher dimensional walks and could result 
in yet another difference between the classical and the quantum walks. We have studied 
the influence of the entanglement and the indistinguishability of the particles on the meeting 
probability. The influence is particularly visible for fermions and in the case of distinguishable 
particles for the case of initial entangled singlet state. Although the meeting probability 
decays faster in these cases the overall meeting probability will still converge to one, as the 
decay is never faster than the threshold |. 
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Conclusions 



Quantum walks are a specialized field on the border between quantum information theory 
and statistical physics which attracted a lot of interest in recent years. A number of novel 
effects have been found and are still under investigation. In the present thesis we contributed 
to these investigations. 

In particular, we extended the concept of recurrence and Polya number to quantum 
walks. The particular measurement scheme employed in our definition preserves the effect of 
the additional degrees of freedom offered by quantum mechanics on the Polya number. We 
developed the tools needed for the analysis of the recurrence nature of quantum walks. The 
actual analysis revealed that quantum walks can be operated in physically different regimes. 
These regimes cover localization as well as ballistic spreading of the walker's wave packets. 
We found that the free parameters we have at hand in a coined quantum walk have a crucial 
impact on its dynamics and are capable of changing its behaviour from recurrent to transient. 
Striking diversity of quantum walks in contrast to classical random walks was pointed out. 
The present results prove the usefulness of the Polya number concept for quantum walks and 
support our expectation of its applicability in related domains. 

The recurrence of quantum walks under the effect of bias was analyzed. For classical 
random walks breaking the symmetry results in immediate turnover from recurrence to tran- 
sience. However, the ballistic nature of the quantum walk is able to compensate for the bias 
and the recurrence can be preserved. We identified the range of parameters for which the 
recurrence behaviour of biased quantum walks on a line diverse from classical random walks. 

Finally, we considered quantum walks with two particles. This makes the additional 
properties offered by quantum mechanics like entanglement or indistinguishability accessible. 
We analyzed the effect of these non-classical features on the meeting probability and pointed 
out the difference from the classical random walk. 

The presented results provide a step in the classification of coined quantum walks, in 
particular on higher-dimensional lattices. We have identified several extreme modes of the 
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dynamics of quantum walks. Our next goal is to exploit the free parameters of the coin 
operator which will allow us to shed light on the connection between these different regimes. 

Within our definition the recurrence of a quantum walk describes the revival of a par- 
ticular quantity, namely the probability at the origin, rather than the revival of a quantum 
state. Nevertheless and quite surprisingly, even full revivals are possible in quantum walk 
settings. This effect is closely related to localization. Indeed, for localizing quantum walks 
the propagator has a non-empty point spectrum which allows for stationary and oscillating 
states. However, the point spectra of the presently known localizing quantum walks are 
rather simple leading only to oscillations with a period of two steps. Finding quantum walks 
with a broader point spectrum will lead to novel features including full and fractional revival 
dynamics. 

Our definition of the Polya number of a quantum walk is connected to a specific mea- 
surement scheme. Needless to say, we can consider schemes where the measurements are 
performed in a different manner and define the Polya number accordingly. It is interesting 
to analyze the influence of various measurement schemes on the recurrence nature of the 
quantum walk. Preliminary results indicate that our definition gives an upper limit for the 
Polya number. 

The meeting problem for two quantum walkers which we have studied presents a step 
towards quantum walks involving many particles. It is certainly worth to analyze other 
various quantities available in multi-particle settings, e.g. the angular correlations among 
the particle's positions. Moreover, up to date the particles performing quantum walk were 
considered non-interacting. To investigate the effect of interactions between the particles on 
the dynamics of quantum walks is one of our next goals. 
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Appendix A 

Recurrence of Random Walks 



In this appendix we review the main results on the recurrence in classical random walks. 
First, we show how the recurrence is related to the probability at the origin. Then we discuss 
the recurrence of unbiased random walks on (^-dimensional lattices. Finally, we analyze the 
recurrence of biased random walks on a line. For a more comprehensive reviews we refer to 
the literature [761175]. 

We begin with the problem analyzed by Polya in 1921 [73j. Consider a particle performing 
a random walk on an infinite ci-dimensional lattice. The particle is initially localized at the 
origin of the lattice. The probability P that the particle returns to the origin during the time 
evolution is called the Polya number of the walk. Random walks are classified as recurrent 
or transient depending on whether their Polya number equals to one, or is less than one, 
respectively. If the random walk is recurrent the particle returns to the origin with certainty. 
On the hand, for transient random walks there is a non-zero probability that the particle 
never returns to its starting point. In other words, there is a non- vanishing probability of 
escape. 

The Polya number of a classical random walk can be defined in the following way |78j . 
Let qo{t) be the probability that the particle returns to the origin for the first time after t 
steps. Since these events are mutually exclusive we can add up their probabilities and the 
series 

oo 

P = (A.l) 

t=l 

gives the probability that at least once the particle has returned to the origin, i.e. the Polya 
number. However, the definition OA.ip is not very practical for determining the recurrence 
nature of a random walk. We can express the Polya number in terms of the probability po(t) 
that the particle can be found at the origin at any given time instant t. Indeed, it is easy to 
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see that the probability at the origin po(t) and the first return probability q (t) fulfills the 
following relations 



Po(0) = 1 

Po(l) = ?o(l) 

Po(2) = <Zo(2)+g (l)Po(l) 

Po(3) = g (3) + g (2)po(l) + go(l)po(2) 

Po(n) = ?o(n) + ?o(n-l)po(l) + --- + ?o(l)Po(n-l)- 

Simply adding all of these equations together might lead to a divergent series. Therefore, we 
first multiply the n-th equation by z n with \z\ < 1. Adding these modified equations we find 
the relation 

F(z) = l + F(z)G(z), (A.2) 
where we have defined the following functions 

oo 

F(z) = £> (n)z" 

n=0 

oo 



G(z) = 5>(n)* B . 



71=1 



Both series are convergent for \z\ < 1. Moreover, the Polya number P can be evaluated from 
the function G(z) by taking the limit z — > 1~ 



P = lim G(^) = Vg 



oo 

,..'")• 

n=l 



From the relation (1A.2j) we express the function Cr(z) in the form 

G(s) = ' ' 



F(s) 

Finally, we take the limit 2 — > 1~ and find the formula 



t=0 

which expresses the Polya number P in terms of the probability at the origin po(t). 
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The recurrence behaviour of a random walk is determined solely by the infinite sum 

oo 

5 = 5>,(f). (A.3) 

t=o 

We find that P equals unity if and only if the series S diverges [78] . In such a case the walk 
is recurrent. On the other hand, if the series S is convergent, the Polya number P is strictly 
less than unity and the walk is transient. The convergence of the series S is determined 
by the asymptotic behaviour of the probability at the origin po(t). Indeed, we find that if 
Po(t) decays faster than t~ l the sum is finite, while if the decay of po(t) is slower the sum 
is divergent. Hence we find the following criterion for recurrence of random walks — the 
random walk is recurrent if and only if the probability at the origin decays like t~ x or slower 
as t approaches infinity. 

In the following we use the above mentioned criterion to analyze the recurrence of biased 
and unbiased random walks. 



A.l Unbiased random walks on 



Let us begin with the unbiased random walk on a line. At each time step the particle has two 
possibilities — it can move to the right or to the left by a unit distance with equal probability 
1/2. The probability distribution generated by such a random walk is easily found to be 

The probability at the origin is thus given by ( for even number of steps 2t ) 

1 (2V 

Using the Stirling's formula 

n\ « V2im y-j (A.4) 
we find that the asymptotical behaviour of the probability at the origin is determined by 

1 



Po(t) 



TYt 



Hence, the series S defined in (1A.3j) is divergent. Consequently, we find that the unbiased 



random walk on a line is recurrent. 
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Recurrence of unbiased random walks on higher- dimensional lattices can be analyzed in a 
similar way [78J . One finds that the asymptotics of the probability at the origin is determined 
by the dimension of the lattice d in the following form 



Po(t) ~ t' 



It follows that the series S determining the recurrence of a random walk is divergent only for 
the dimensions d = 1, 2 and convergent for d > 3. We conclude that the random walks on a 
line and in the plane are recurrent while higher-dimensional random walks are transient, the 
result originally found by Polya in 1921 [73] . 

Concerning the value of the Polya number for the transient case Montroll [H] showed 
that for the dimensions d > 2 the following relation holds 

1 



P(d) = 1 - 



u(d) 



where u(d) can be expressed in terms of an integral of the modified Bessel function of the 
first kind 



u(d) 



o 



e' l dt. 



However, the closed form of the function u(d) is known only for d = 3 due to the Watson's 
triple integral [82] with the result 

-(3) = (L) r (*) r (L) r (") „ L51 , 



32vr 3 V 24 / V 24 / V 24 / V 24 / 

For higher dimensions d > 3 one has to evaluate the integral numerically. We present an 
overview of the numerical values of the Polya number [71] for a different dimensions d in 
Table [Ql 



A. 2 Biased random walks on a line 

Let us consider biased random walks on a line. The bias can be introduced in two ways - 
the step in one direction is greater than in the other one and the probability of the step to 
the right is different from the probability of the step to the left (see Figure POT) . 

Consider a random walk on a line such that the particle can make a jump of length r to 
the right with probability p or make a unit size step to the left with probability 1 — p. As we 
have discussed a random walk is recurrent if and only if the probability to find the particle 
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Dimension 


Polya number 


3 


0.340537 


4 


0.193206 


5 


0.135178 


6 


0.104715 



Table A.l: Polya number of a random walk on Z d in dependence of the dimension d. 



1 - p V 




-2-10 1 r r + 1 

Figure A.l: Schematics of the biased random walk on a line. The particle can move to the 
right by a distance r with the probability p. The length of the step to the left is unity and 
the probability of this step is 1 — p. 



at the origin at any time instant t does not decays faster than t 1 . This probability is easily 
found to be 



P0[t) = (1 — p)r+lpr + l 



t 

tr 
r+1 



With the help of the Stirling's formula ( 1A.4j) we find the asymptotical behaviour of the 
probability at the origin 



Po(t) 



r + 1 
V2nrt 



'1 \ — 2i_ _j_r + 1 

K l — P) r+lpr+l __ 

f r+1 



-I t 



The asymptotics of the probability po(t) therefore depends on the value of 

M \ — — ^r + 1 

q = (1 — p)r + lpr+l r . 

nr r+1 

Since q < 1 the probability po(t) decays exponentially unless the inequality is saturated. 
Hence, the random walk is recurrent if and only if q equals unity. This condition is satisfied 
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for 




(A.5) 
gth of 




the step. 



This result can be well understood from a different point of view, as we illustrate in 
Figure IA.2I The spreading of the probability distribution is diffusive, i.e. a ~ y/i. The 



neighborhood the probability decays exponentially. Therefore for the random walk to be 
recurrent the origin must lie in this a neighborhood for all times t. However, if the random 
walk is biased the mean value of the position (x) varies linearly in time, thus it is a faster 
process than the spreading of the probability distribution. In such a case the origin would lie 
outside the a neighborhood of the mean value after a finite number of steps leading to the 
exponential asymptotic decay of the probability at the origin po(t). Hence, the random walk 
is recurrent if and only if the mean value of the position equals zero. Since the individual 
steps are independent of each other the mean value after t steps is simply a t multiple of the 
mean value after single step, i.e. 



probability in the a neighborhood of the mean value (x) behaves like t 2 while outside this 



(x(t)) =t(x{l)) =t\p(r + l) -1]. 



We find that the mean value equals zero if and only if the condition (1A.5j) holds. 
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Figure A. 2: Spreading of the probability distribution versus the motion of the mean value 
of a biased classical random walk on a line. While the spreading is diffusive (er ~ y/i) the 
mean value propagates with a constant velocity ((x) ~ t). The probability inside the a 
neighborhood of the mean value (x) behaves like t". On the other hand, outside the a 
neighborhood the decay is exponential. Hence, if the mean value (x) does not vanish the 
origin of the walk leaves the a neighborhood of the mean value. In such a case the probability 
at the origin decays exponentially and the walk is transient. 
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Appendix B 

Recurrence Criterion for Quantum 
Walks 

Let us prove that the recurrence criterion for quantum walks is the same as for random walks, 
i.e. the Polya number equals one if and only if the series 

oo 
t=0 

diverges. 

According to the definition of the Polya number Eq. (11. 5p for quantum walks we have to 
prove the equivalence 

+00 

P = 11 (1 - Po {t)) = ^ S = +00. 
t=i 

We note that the convergence of both the sum S and the product P is unaffected if we omit 
a finite number of terms. 

Let us first consider the case when the sequence Po(t) converges to a non-zero value 
< a < 1. Obviously, in such a case the series S is divergent. Since po(t) converges to a we 
can find for any e > some to such that for all t > to the inequalities 

1 — a — e < 1— Po(t) < 1 — a + e. 

hold. Hence, we can bound the infinite product 

lim (l-o-e)*<P< lim (l-o + e)*. (B.l) 

t— >+oo t— >+oo 

Since we can choose e such that 

|1 — a ± e\ < 1, 
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we find that limits both on the left-hand side and the right-hand side of Eq. (IB . 1 [) equals 
zero. Hence, the product P vanishes. 

We turn to the case when po(t) converges to zero. We denote the partial product 

n 

Pn = l[(l-P0(t)). 

t=l 

Since 1 — Po{t) > for all t > 1 we can consider the logarithm 

lnP n = ^ln(i-j9o(t)) (B.2) 

i=l 

and rewrite the infinite product as a limit 

P = lim e lnTn . (B.3) 



Since po(t) converges to zero we can find some to such that for all t > to the value of po(t) is 
less or equal than 1/2. With the help of the inequality 

— 2x < In (1 — x) < —x 

valid for x G [0, 1/2] we find the following bounds 

n n 

-2j2Po(t) < lnP n <-J2po{t). 

t=\ t=i 

Hence, if the series S is divergent the limit of the sequence (in P n ) is — oo and according to 
Eq. ( 1B.3|) the product P vanishes. If, on the other hand, the series S converges the sequence 

(in P n ) n _ 1 is bounded. According to Eq. (IB.2j) the partial sums of the series ^ In (1 — po(t)) 

n ~ t=i 
are bounded and since it is a series with strictly negative terms it converges to some negative 

value b < 0. Consequently, the sequence (lnP n )^ =1 converges to b and according to Eq. 

(IB.3P the product equals 

P = e b > 0. 

This completes our proof. 
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Appendix C 

Method of Stationary Phase 



In order to determine the recurrence nature of a quantum walk one has to analyze the 
asymptotic behaviour of the probability at the origin. As we have shown in Section 11.41 the 
probability amplitude of the particle being at the origin of the quantum walk after t steps is 
given by a sum of integrals of the form 



The recurrence of a quantum walk is determined by the asymptotics of such integrals. The 
method of stationary phase is a suitable tool for such analysis. 

In the following we briefly review the main concepts of the method of stationary phase. 
First, we treat the one- dimensional integrals. Then we turn to the multivariate integrals. We 
find that the crucial contribution to the integral ( IC.ip as t approaches infinity arises from the 
stationary points, i.e. the points where the derivative of the phase w(k) vanishes. We discuss 
how the amount of stationary points and the "flatness" of the phase at the stationary point 
influences the asymptotic behaviour of the integral I(t). For a more comprehensive analysis 
we refer to the literature [531 ED] • 

C.l One-dimensional integrals 

Let us begin with the one- dimensional integral of the form 

b 



where / and u are smooth functions and u is real-valued. We see that in the region of k 
where u(k) changes considerably the exponential e % oscillates rapidly as t approaches 




(C.l) 



v 





a 
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infinity. Assuming that the function / is slowly varying compared to these rapid oscillations 
we find that this region of integration does not contribute significantly to the integral I(t). 
Obviously, the most important contributions to the integral (IC.2j) arise from the regions where 
the oscillations of the exponential are least rapid, which occur precisely at the stationary 
points k of the phase u 

i / , % (ho 

*/%) = Tk = 0. 

feo 

The "flatness" of the phase at the stationary point determines the order of this contribution — 
the more derivatives of the phase vanishes at the stationary point the slower the contribution 
decays as t approaches infinity. Here we assume that the function / is non-zero at the 
stationary point, otherwise the contribution to the integral I(t) vanishes. 

C.l.l No stationary points 

Let us first consider the case when the phase u has no stationary points inside the integration 
domain. Then there exists e > such that 

\u'(k)\ >e 

for all k. Performing the integration in (1C.2j) per parts with 



U (k) = v'(k) = %tj{ky w < fc >*, 

u'(k) - ^f'(kW(k)-f(kW>(k) v(k) _ eiuiuyt 

u[k) - it u'(ky ' v{k) ~ e 



we find that I(t) can be expressed in the form 



7 ^ - It 



f(k) 
u'(k) 



b 



1 f r j „( k) t f'(kW(k)-f(kW'(k) 

It J 6 Z7W { } 



a 



We see that I(t) decays at least like t -1 as t approaches infinity. Moreover, the second term 
in (1C.3j) has the same form as the original integral (1C.2j) . Hence, if in addition the first term 
in (IC.3j) vanishes, e.g. if the function / equals zero at the boundaries of the integration 
domain, we find by repeated integration per parts that I(t) decays faster than any inverse 
polynomial in t. 
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C.1.2 First-order stationary points 

We turn to the case of u having a single stationary point coinciding with the left endpoint 
of the interval ko = a ( any integral where the phase has more than one stationary point can 
be decomposed into a sum of such integrals ) and assume that the stationary point is of the 
first order, i.e. u'(a) = but u"(a) ^ 0. We then expand the phase into a Taylor series 

u(k) ~u(a) + ^^-(k-a) 2 

around the stationary point ko = a. Since we assume that the function / is slowly varying 
we put it equal to its value at the stationary point f(k) m /(a). With these estimations we 
find 

b 

I{t) ~ /(ay " (a) * J e^^-^dk. (C.4) 

a 

Let us estimate the remaining integral. We substitute for y = k—a and extend the integration 
domain to [0, +oo) 

b b—a +oo 

f^-^ dk= [ e ^y 2t dy* f e^ H dy. 



a 

This is the familiar Fresnel integral 



F( 7 )= / e^ y2 dy= Ly.| 7 |-| e »ign7f. (c.5) 







Inserting this result into the expression ( 1C.4j) we finally arrive at the formula 



I(t) 



t~2 (C.6) 



K(a)|. 

describing the behaviour of the integral 1(f) for large values of t. The plus (minus) sign in 
the exponential in (1C.6P corresponds to the second derivative of the phase at the stationary 
point u"(a) being positive (negative). To conclude, we find that if the phase u{k) has a 
stationary point of the first order the integral lit) decays like t _1//2 as t approaches infinity. 

C.1.3 Higher-order stationary points 

We close this section by the analysis of the integral I(t) when the phase uo has a stationary 
point k = a of the order of p — 1, i.e. 

io'(a) = u"{a) = ... = uj^ia) = 0, oo^(a) ^ 0. 
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In such a case the Taylor expansion of the phase reads 

u(k)~u(a) + ^-^(k-a) p . 
pi 

Performing similar approximations as above we find 



b 

I(t) ~ /(a)^ w <°>* [ j^^-^dk. (C.7) 



In the remaining integral we substitute for y = k — a and extend the upper limit of the 
integration to +00 

b b—a +00 

e p ! k ' ak = e p- ay w / e ?> ! ay. 
a 00 

We find a generalization of the Fresnel integral (1C.5P which is readily evaluated 

+00 



F p ( 7 )= / e*^dy= -^^\ 7 \-i e imsn ^. (CJ 







P 

Finally, inserting this result into the Eq. (10. T|) we arrive at the estimation 



I(t) 



/W P VM p) (a)l 



t~p, (C.9) 



where the plus (minus) sign corresponds to positive (negative) value of wW(o). From ffCTQ^ 
we find that the contribution of the stationary point of the order p — 1 to the integral I(t) 
behaves like t~ l l p as t approaches infinity. The flatness of the phase at the stationary point 
reduces the rate at which the integral I(t) decays. 

C.2 Multivariate integrals 

We turn to the asymptotic analysis of the multidimensional integrals of the form 

I(t) = [ e l " (k)i /(k)dk. (CIO) 



v 
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We assume that both functions w(k) and /(k) are smooth and u is real- valued. Similarly 
to the one-dimensional case, the main contribution to the integral arise from the stationary 
points of the phase u, i.e. points k where the gradient of uj vanishes 



Vu;(k)| k=ko = 0. 

As in the previous Section we approximate the phase around the stationary point by the 
Taylor expansion. In addition, we have to change the coordinates in such a way that the 
resulting integral factorizes into a product of one- dimensional integrals. Each of the 1-D 
integrals can be estimated by means provided in the previous Section. 

In the following we review the main results of the asymptotics of (IC.lOp in dependence 
on the properties of the phase u(k). For a more detailed analysis we refer to the literature 
[831 [80]. 

C.2.1 No stationary points 

Let us begin with the case when the gradient of u is non-vanishing inside the integration 
domain V. From the divergence theorem we find 

I(t) = -- J{u • n)e i ^ds + l - J (V ■ uy ut dk, (C.ll) 

dV V 

where dV is the boundary of V, n is the unit vector normal to the boundary and the vector 
function u(k) is given by 

The expression (IC.llj) indicates that I(t) decays at least like t™ 1 . Suppose that the function 
/(k) vanishes smoothly on the boundary of V. In such a case the contour integral in (IC.llj) 



equals zero. The remaining volume integral in (1CJ.11)> is of the same kind as the original 
integral I(t). Hence, by repeating the same procedure as above we find that the integral I(t) 
decays faster than any inverse polynomial in t. 

C.2.2 Non-degenerate stationary points 

We turn to the case when the phase w(k) has a single stationary point ko inside the integration 
domain. We assume that k is non-degenerate, i.e. the Hessian matrix evaluated at the 
stationary point 



i^(ko) 



d 2 cu 
dkidkj 



(C.12) 

k=k 
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is regular. We expand the phase around the stationary point into the second order 
w(k) ~ w(k ) + - ( k i ~ k oi) H i,j( k o) ( k j ~ k oj) ■ 



Assuming that /(k) is slowly varying we can evaluate it at the stationary point and extract 
it from the integral. Substituting for 

k, = k — ko 

and extending the integration from the finite volume V to R n we arrive at the following 



estimation of the integral fIC.lOj) 



lit) ~ /(k )e J / exp KiH^Kjt dn. 



(C.13) 



The integral in (IC13P can be reduced into the product of n one-dimensional Fresnel integrals 
(lC.5p . Indeed, the Hessian matrix (10. 12[) is real and symmetric since we assumed oo(k) to be 
smooth. Hence, it can be diagonalized with the help of the orthogonal matrix O. In the new 
coordinate system 

^ = ^0^^ (C.14) 



the bilinear form in (IC.13j) is given by the sum of purely quadratic terms 

(ko)/^- = ^Ai(k )/4 



where Aj(k ) are eigenvalues of the Hessian matrix (IC.12P at the stationary point k . Since 
the matrix O is orthogonal the change of coordinates (10.14[) has a unit Jacobian. Hence, 
using the substitution (1C.14j) we decompose the integral in (1C.13I) into the product of one- 
dimensional Fresnel integrals 



/ 

JR." 



exp ( o H ' l, 'j( k '>)'< J l J dn = Jj£ / exp ( -Aj^ko)^ ) dfij, 



3=1 



JR. 



which are readily evaluated with the help of (IC.5I) . Finally, we arrive at the following ap- 
proximation of the integral (IC.lOp 



I(t) 



/(k )e 



\det Hi j(k ) 



t"2 



(C.15) 



where ^(k ) is the sum of the signs of the eigenvalues of the Hessian matrix 



K k o) = 5^ signA j( k °)- 

3 

We find that contribution from the non-degenerate stationary points to the n-dimensional 
integral ( IC.lOp is of the order of t~ n / 2 . 



C.2.3 Continuum of stationary points 



We close this Appendix by briefly discussing the asymptotic scaling of the integral ( IC.lOj) 
when the phase oo(k) has a curve of stationary points 7, i.e. 

Vk e 7 Vw(k) = 0. 



Without loss of generality we assume that w(k) = at the stationary curve 7 which is 
considered to be smooth and without any loops. Moreover, we restrict ourselves to two- 
dimensional integrals, i.e. n = 2. As shown in [80J, Chapter VIII. 9, the main contribution of 
the continuum of stationary points to the asymptotic expansion of the integral f lC.lOj) is 



1(f) 



27re 4 4 



f(h(s),k 2 (s)) 



d 2 ui I d' 2 u. 
dk'f dk^ 



ds 



where s is the parametrization of the curve 7. We find that in comparison with the case 
of the isolated non-degenerate stationary point analyzed in Section IC.2.21 the continuum of 
stationary points has reduced the decay of the integral I(t) by a factor of square-root. 
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Appendix D 
Meeting Problem 



In this Appendix we analyze the meeting problem in classical and quantum walk. We derive 
analytical formulas for the asymptotic behaviour of the meeting probability. 

D.l Meeting problem in the classical random walk 

Let us define the meeting problem on the classical level. We assume two particles which 
in each step of the process can perform randomly a step to the left or to the right on a 
one dimensional lattice labelled by integers. Initial distance between the two particles is 
2d, because for odd initial distance the two particles never meet, due to the transitional 
invariance we can assume that one particle starts from the origin and the other one in the 
vertex 2d. We assume complete randomness, i.e. the probabilities for the step right or left are 
equal. We ask for the probability that the two particles meet again after t steps either at a 
certain position m or we might ask for the total probability to meet (the sum of probabilities 
at all of the possible positions). A simple analysis reveals that the probability to meet at a 
certain position m equals 




The total probability that the two particles are reunited after t steps reads 




which simplifies to 




(D.l) 
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To obtain the asymptotic behavior of the meeting probability we approximate the single 
particle probability distribution by a gaussian 

(x -2d) 2 " 



P d (x,t,d) 



1 



exp 



'id ' V 2t 
which leads to the following estimate on the meeting probability 



+00 



f 1 / d 2 

M cl (t, d) « J P d (x, t, 0)P d (x, t, d)dx = -j= exp 1-— 



Finally, for a fixed initial distance d we get the long-time approximation for t > d 2 



D.2 Meeting problem in the quantum walk 

Let us derive analytical formulas for the meeting probabilities in the quantum case. We 
consider the following initial states: 

(i) right for the first particle and left for the second 

^x(0,2rf,0) = 1, 

(ii) symmetric initial conditions l/y/2(\L) + i\R)) for both 



/ 1 \ 



#),2d,0) = ~ 



V- 1 / 

(Hi) left for the first particle and right for the second 

*P LR (0,2d,0) = l. 

For t > \/2d we consider the slowly varying part of the single particle probability distri- 
bution derived in [51] which has the form 

2 



slow 



x,t) 



2x 2 
t 2 
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if the initial coin state was \L) or \R), while for the symmetric initial condition it reads 



p(S) 

slow 



2x 2 
t 2 



^(i-f)V 1 - 

We then estimate the sums in ( I4.19|) defining the meeting probabilities by integrals 



M RL (t,d) 



t 

2 f dx 



- h \j (i-f)(i+^) A /r^s A /i-2^ 



M s {t,d) 



2d--!- V x (A iT t >V t 2 V t2 

2 f dx 



rt \4 d - ff )d - ^ty 1 - 2 'H l - 2 '~- 



V2 

2 f dr 
M LR (t,d) « -= / = (D.2) 

which can be evaluated in terms of elliptic integrals. Notice that the integrals diverge for 
d = 0, i.e. for the case when the two particles start at the same point. For now we suppose 
that d > 0. The formulas (ID.2j) can be expressed in the form 



M RL (t, d) « F + |2(t - d)(t - (4 - 2 v / 2)d)^(a) + 

+V2 ((t - (4 + 2v / 2)d)(t - (4 - 2 v / 2)rf)n(6+|a) - t 2 n(c + |a) ) } 



M s (t,d) « |l6rf(t 2 -rf 2 )(t + (4 + 2v / 2)d)(t- (4 - 2V2)d)A"(a) + 

+v / 2(t + (4 + 2V2)tf)(* - (4 + 2\/2)d)(t + (4 - 2^2)^) x 
x(t - (4 - 2\/2)d) f(t + d)n(6+|a) + (t- d)n(6_|a) ) - 

-v 7 ^ 2 ((t + d){t + (4 + 2v / 2)d)(t + (4 - 2v / 2)rf)n(c+|a) + 
+(t - d){t - (4 + 2V2)d)(t - (4 - 2v / 2)d)n(c_|a) ) | 

M LR (t,d) « F_ {2(t + d)(t + (4 + 2^)^)^(0) - 

-v 7 ^ + (4 + 2v / 2)d)(t + (4 - 2 v / 2)c/)n(6_|a) -t 2 n(c_|a) 



(D.3) 
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Here K(a) is the complete elliptic integral of the first kind and n(x|a), n(y|a) are the complete 
elliptic integrals of the third kind (see e.g. [SI], chapter 17). The coefficients a, b±,c± and 
F± are given by 



2t 



7r 2 d(t T d)(t(2 + y/2) =F Ad)(t(2 - y/2) T 4c/)) 
2d 2 

{l±y/2){t-y/2d) 

d(V2~T2) 
(t(y/2^2) + 4d)(t-y/2d) 
y/2d(t{y/2±2) -Ad) 



Let us analyze the asymptotic behavior of the meeting probability. We begin with the 
observation that the coefficients at the highest power of t with the elliptic integrals of the 
third kind are the same but with the opposite signs for n(6|a) and IT(c|a). Moreover, b± 
and c± goes like —t as t approaches infinity, and thus all of the II functions have the same 
asymptotic behavior. Due to the opposite sign for n(6|a) and II(c|a) the leading order terms 
cancel and the contribution from this part to the meeting probability is of higher order of 
\jt compared to the contribution from the complete elliptic integral of the first kind K{a). 
The asymptotic of the function K(a) is given by 

rfv /21n 
K(a) « 

Inserting this approximation into (ID. 3}) we find that the leading order term of the meeting 



probability in all the three studied situations is given by 

In 

M D (t,d) 
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Part II 

Factorization with Exponential Sums 
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Introduction 



Factorization of integers is a famous NP problem jSU [85] and the difficulty to decompose a 
number into prime factors lies at the heart of several encryption schemes [861 EZ]- However, 
Peter Shor found [88] that a quantum computer is capable of finding factors of a given 
number efficiently The fundamental advantage of the Shor's algorithm compared to the 
classical algorithms is the massive use of quantum parallelism and entanglement. On the 
other hand, the physical realizations of the Shor's algorithm are very challenging and are so 
far limited to a proof of principle experiment [89J. 

Recently, several schemes for integer factorization based on Gauss sums [90, ED E2] were 
proposed [931 EH ESI ESI EH M, EH1 HDD]. For a review see e.g. [M]. In contrast to the 
Shor's algorithm, factorization using Gauss sums consists of a feasible factor test based on 
a classical interference scheme. The proposals employ the so-called normalized truncated 
Gauss sum 



Here N is the number to be factored, £ is a trial factor and M is the truncation parameter. 
The capability of Gauss sums to factor numbers stem from the fact that the signal - the 
absolute value of the Gauss sum which is measured in the experiment, attains the maximal 
value only for a factor. For non-factors destructive interference yields a small signal. In the 
most elementary approach we have to perform this factor test for every number smaller than 
"v/iV. As a consequence the method scales as yN and is therefore exponential. On the other 
hand, the physical realizations of Gauss sums are less demanding than the implementations of 
the Shor's algorithm. Indeed, recent experiments based on NMR |102[ 11031 HE], cold atoms 
|105] , ultra-short pulses |1061 I107j [I08J and Bose-Einstein condensate [109] have successfully 
demonstrated the possibility to find the prime factors of up to 17-digit numbers. 

In the NMR settings [1021 11031 I104j a sequence of RF pulses with linearly increasing 
relative phase shifts is applied to the ensemble of nuclear spins. After each pulse the echo, 
i.e. the polarization of the spins, is measured. Finally, all echoes are summed and for the 
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proper choice of relative phase shifts of the RF pulses the resulting signal has the form of 
the Gauss sum. 

The experiment with cold atoms presented in [105] employs two long-living hyperfine 
ground states of rubidium. The atoms are launched by a system of magneto-optical traps and 
prepared in the atomic ground state by appropriate pulse sequence. After the preparation the 
atoms interact with a sequence of Raman pulses driving a transition between the hyperfine 
states. Similarly to the NMR experiments the individual pulses have to be properly phase 
shifted. Finally, after all pulses are applied a fluorescence detection measures the populations 
in both hyperfine states. The sum of these interference signals determines the Gauss sum. 

In [106^ 1107} 1108] the Gauss sum is implemented by a sequence of shaped femtosecond 
laser pulses. Individual laser pulses are properly phase shifted by a complex spectral mask. 
The interference produced by the pulse train is analyzed with a spectrometer. Due to the 
temporal Talbot effect the frequency component of the electric field is determined by a Gauss 
sum. 

The experiment [109] uses diffraction of the BEC on an optical lattice. One of the beams 
which creates the optical lattice is designed with a specific phase jumps. The pulse separates 
the atoms in the BEC into different momentum orders. In the absorption image a diffraction 
pattern determined by the Gauss sum is observed in which high-momentum atoms represent 
factors and low-momentum atoms represent non-factors. 

As we have mentioned, the signal for a non-factor is suppressed and its value depends 
on the number of terms in the Gauss sum. In the experiment we have to take into account 
the limited resolution of the measured signal. Hence, to be able to distinguish factors from 
non-factors we have to add a sufficient number of terms in the Gauss sum. However, in 
all experiments performed so-far the individual contributions to the Gauss sums are created 
by individual pulses. Hence, the total number of terms in the Gauss sum is limited by the 
decoherence time of the system used in the experiment. Because of these two antagonistic 
effects we have to find conditions under which the algorithm based on Gauss sums successfully 
finds the factors of a given number N. We answer these questions in the following Chapters. 

Chapter [5] deals with truncated Gauss sums and is based on [VI] . We find that the 
truncated Gauss sums offer good discrimination of factors from non-factors since the gap 
between their corresponding signals can reach a value of almost 30%. Moreover, we show 
that to reach such a gap the number of terms in the Gauss sum M we have to add, i.e. 
the number of laser pulses we have to apply in the experiment, has to be of the order of the 
fourth-root of N. The total number of the resources needed for the success of the factorization 
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scheme based on the truncated Gauss sum is thus 



n 



'N-VN 



In Chapter E] which is based on |VIIj we extend the idea of factorization of integers from 
Gauss sums to exponential sums of the form 



M 



m=0 



2ni m? - 



Here the power of the phase is no longer quadratic like in the case of the Gauss sum but 
is given by a positive integer j. The faster growth of the phase results in the reduction of 
the number of terms M that has to be added to the 2j-th root of N. The total number of 
resources necessary to factorize a number N using exponential sum is given by 



3 + 1 



Hence, we can save experimental resources by applying exponential sums with larger value of 
j. On the other hand, the gap between the signals of factors and non-factors shrinks as the 
power of the phase j increases. This can make the experimental data inconclusive, unless a 
sufficient resolution is guaranteed. We summarize our results in the Conclusions. 
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Chapter 5 

Factorization with Gauss sums 



Introduction 

Gauss sums [90j [9TJ [92] play an important role in many phenomena of physics ranging from 
the Talbot effect of classical optics [110] via the curlicues emerging in the context of the semi- 
classical limit of quantum mechanics |111[ 1112] , fractional revivals |1131 1114] and quantum 
carpets |115j to Josephson junctions |116j . Moreover, they build a bridge to number theory, 
especially to the topic of factorization. Indeed, they can be viewed as a discrimination 
function of factors versus non-factors for a given natural number. The essential tool of this 
factorization scheme jHH] is the periodicity of the Gauss sum. 

Usually Gauss sums extend over some period which leads to the complete Gauss sum. 
However, recent experiments based on NMR |102[ I103[ 1104] . cold atoms |105j . ultra-short 
pulses |106[ I107[ 1108] and Bose-Einstein condensate [109] have demonstrated the possibility 
of factoring numbers using a truncated Gauss sum, where the number of terms in the sum 
is much smaller than the period. Therefore, factorization with truncated Gauss sums offers 
enormous experimental advantages since the number of terms is limited by the decoherence 
time of the system. In the present Chapter we address the dependence of the number of 
terms needed in order to factor a given number. In particular, we find an optimal number 
of terms which preserves the discrimination property and at the same time minimizes the 
number of terms in the sum. 

In order to factor a number N we analyze the signal, i.e. the absolute value of the 
Gauss sum, for integer arguments I = 1, . . . , |_V^Vj- We call the graphical representation 
of the signal data factorization interference pattern. In order to gain information about the 
factors of N we analyze the factorization interference pattern: Whenever the argument I 
corresponds to a factor of N we observe the maximal signal value of unity. For most non- 
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factor arguments this signal value is significantly below unity. However, for ghost factors we 
observe signal values close to unity even though these arguments do not correspond to an 
actual factor of N. Thus ghost factors spoil the discrimination of factors from non-factors in 
such a factorization interference pattern. Fortunately, ghost factors can be suppressed below 
a given threshold by extending the upper limit of the summation in the Gauss sum. This goal 
of completely suppressing all ghost factors provides us with an upper bound on the truncation 
parameter. This upper bound represent a sufficient condition for the success of our Gauss 
sum factorization scheme. The analysis of the number of ghost factors evaluated by the ghost 
factor counting function g(N, M), which depends on the number to be factorized N and the 
truncation parameter M, reveals that this upper bound on the truncation parameter is also 
a necessary condition for the success of our Gauss sum factorization scheme. 

The Chapter is organized as follows: We first briefly review in Section [57X1 the central idea 
of the factorization scheme based on the Gauss sums. In particular, we introduce complete 
and truncated Gauss sums and compare the resources necessary to factor a given number 
N. We find the first traces of ghost factors in the factorization interference pattern based on 
the truncated Gauss sum. Since the truncation of the Gauss sum weakens the discrimination 
of the factors from non-factors, we dedicate Section 15.21 to deriving a deeper understanding 
of this feature. We find four distinct classes of arguments i which result in utterly different 
behaviours of the truncated Gauss sum. Rewriting the truncated Gauss sum in terms of the 
curlicue sum allows us to identify the class of problematic arguments - the ghost factors. 
Moreover, we identify a natural threshold which separates factors from non-factors. For 
a rigorous argument we refer to Appendix [El In Section 15.31 we obtain an upper bound 
on the truncation parameter of the Gauss sum needed to suppress the signal of all ghost 
factors below the natural threshold. Ghost factors appear whenever the ratio of the number 
to be factored and a trial factor is close to an integer. This fact allows us to replace the 
Gauss sum by an appropriate Fresnel integral. From this expression we find the scaling law 
M ~ \fN for the truncation parameter M, which represents the sufficient condition for the 
success of our Gauss sum factorization scheme. We discuss the applicability of the Fresnel 
approximation in the Appendix [F] Finally, we analyze the ghost factor counting function 
in Section 15.41 and show that the fourth-root law is also necessary for the success of our 
factorization scheme, even if we relax the threshold value or allow limited error tolerance. 
We conclude in Section [531 
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5.1 Factorization based on Gauss sums: appearance of 
ghost factors 



To start our analysis we first consider the complete normalized quadratic Gauss sum 



A^ 1 \£) = ]^e W (2 mm ^), (5.i; 

m— n V / 



m=0 

which is frequently used in number theory. Here N is the integer to be factorized and the 
integer argument £ scans through all numbers from 1 to Lv^J f° r factors of N. If £ is a 
factor then all terms in the sum contribute with a value of unity and thus the resulting signal 
value l \£)\ is one. However, for non-factor arguments the signal value is suppressed 



VN 

^£ = -y/N(VN- 1) w -N (5.2) 



considerably as illustrated on the left in Figure 15.11 Thus the absolute value of the Gauss 
sum allows one to discriminate between factors from non-factors. 

Factorization based on the complete Gauss sum ( 15.1 j) has several disadvantages. First of 
all, the limit of the sum depends on the trial factor I. Thus the number of terms in the sum 
increases with £ up to \/N. Hence, to obtain a complete factorization interference pattern in 
total 

m 

£= 

2 ' 2 

i=\ 

terms have to be added. 

In the recent experimental demonstrations pH Q3S1 QUI [HH QUI [Ml CEDE! EE] of our 
Gauss sum factorization scheme the number of terms in the sum translates directly into the 
number of pulses applied onto the system, or the number of interfering light fields. Due 
to the decoherence it is favorable to use as few pulses as possible. Hence the experiments 
employ a constant number M of pulses for each argument £ to be tested. Thus the resulting 
signal is of the form of a truncated Gauss sum 

Af\£) = j^-j J2 e*P ( 2mm2 j) » ( 5 - 3 ) 

771=0 ^ ' 

rather than a complete Gauss sum of (15. ip . Hence we have to add 

J^M = M ■ VN (5.4) 
e=i 

terms to obtain the factorization pattern with the truncated Gauss sum. With this fact in 
mind we treat the number of terms in the Gauss sum as a resource for this factorization 
scheme. 
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The experiments impressively demonstrate that the truncated Gauss sums are also well 
suited to discriminate in the factorization interference pattern between factors from non- 
factors, even though the summation range does not cover a full period. As a drawback we 
find that the signal value at non-factor arguments is not suppressed as well as in the case of 
the complete Gauss sum. 

In order to illustrate the effect of truncating the Gauss sum we compare in Figure 15.11 
the factorization interference patterns for the complete Gauss sum A^~ 1} (e) (left) and for 
the truncated Gauss sum A { x\l) (right). In a first guess we chose the truncation parameter 
M = In N to depend logarithmically on the number to be factorized. It is remarkable that 
the small number M = 16 of terms in the truncated Gauss sum is sufficient to reveal the 
factors of a seven-digit number like N = 9624687. On the other hand we observe a number 
of data-points with signal values close to one (stars), for example at the argument I = 2555. 
In an experiment such points might lead us to wrong conclusions in the interpretation of a 
factorization interference pattern. Thus we call arguments £ resulting in such critical values 
of the signal ghost factors. 

5.2 Classification of trial factors 

The frequency of appearance of ghost factors is the central question of our study. Indeed, 
how many terms in the truncated Gauss sum are needed in order to suppress the occurrence 
of ghost factors. However, we first need to identify the class of arguments which results in 
ghost factors. 

Figure 15.11 already indicates that there are different classes of trial factors: (i) factors 
with constant value of |.A^(£)|, (ii) typical non-factors at which already few terms M are 
sufficient to suppress the value of \A}f considerably, (iii) ghost factors at which a larger 
summation range is needed to suppress the value of |*4.^(£)|, and finally (iv) threshold 
non-factors at which the value of (•Ol cannot be suppressed by increasing M. 

To illustrate this we plot in Figure I5T21 the truncated Gauss sum of (15.31) for iV = 9624687 
and various arguments i characteristic for each one of the class (i-iv) as a function of the 
truncation parameter M. 

To which class of arguments (i-iv) the given I belongs is determined by the relation 
between the argument i and the number we are factorizing N, namely on the value of the 
fraction 2N/£ which enters the Gauss sum (15.31) . Indeed, for the number N = 9624687 and 
the arguments I used in Figure [531 we find the following: (i) for a factor i = 919 the fraction 
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3- 919 





Figure 5.1: Influence of the truncation parameter M of the incomplete Gauss sum (£) 
defined in (I5.3P on the contrast of the factorization interference pattern for the example 
iV = 9624687 = 3 • 919 • 3491. The upper picture shows the corresponding pattern for the 
complete Gauss sum defined in (15.11) where M = £ — 1 . For the lower plot we have truncated 
the Gauss sum after M = InN = 16 terms. At factors of N indicated by vertical lines the 
Gauss sum assumes the value of unity marked by black diamonds. The complete Gauss sum 
enjoys an impressive contrast due to a suppressed signal value at all non- factors. However, 
also the truncated Gauss sum with a relatively small number of terms allows to discriminate 
factors from non-factors. However, we also observe several ghost factors marked by stars. 
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N = 9624687 = 3 • 919 ■ 3491 




Figure 5.2: Emergence of four different classes of arguments £ of the truncated Gauss sum of 
(15. 3p from its dependence on its truncation parameter M for N = 9624687 = 3 -919 -3491. We 
show the signal value (£) | for four arguments i as a function of the truncation parameter 
M. For factors of N, such as £ = 919 depicted by the black diamonds, the signal is constant 
and equal to unity. For typical non-factors, such as £ = 14 depicted by the gray dots, the 
signal is suppressed considerably already for small values of the truncation parameter M. 
However, for ghost factors, such as £ = 2555 depicted by stars, much more terms in the sum 
(I5.3p are needed to suppress the signal. For arguments such as £ = 12 depicted by the black 
triangles, the signal levels off at a non- vanishing threshold determined by 1/ a/2 and it is 
impossible to suppress it by further by increasing the truncation parameter M. 



2N/£ is an even integer, (ii) for a typical non-factor £ = 14 the fraction 2N/£ is close to an 
odd integer, (iii) for a ghost factor £ = 2555 the fraction 2N/£ is close to an even integer, (iv) 
for a threshold non-factor £ = 12 the fraction 2N/£ is an even integer plus one-half. Thus we 
see that the class of £ is given by the fractional part of the fraction 2N/£. Hence, in order to 
bring out these classes most clearly, we represent the truncated Gauss sum (15.31) in a different 
form. For any argument £ we decompose the fraction 2N/£ into the closest even integer 2k 
and the fractional part p(N,£) = p/q with \p\ < 1 and p, q being coprime, i.e. 

p(N,£) = ™-2k. (5.5) 
Since exp (2tii m 2 ■ k) — 1 the Gauss sum (15.31) reads 

Af\£) = s M {p{N4)) (5.6) 
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where we have introduced the normalized curlicue function |111] . [TT2"] 

1 M 

SM ^ " M+ 1 ^ 6XP ( i7im2r ) ( 5 - 7 ) 

m=0 

which we consider for a real argument r with — 1 < r < 1. 

The connection (15. 6p between the truncated Gauss sum A { x \l) defined in and the 
normalized curlicue sum sm{t) for a given N is established by the fractional part p(N,£) of 
the fraction 2N/£. Indeed, factors of N correspond to p = 0. All other values of p correspond 
to non-factors. In particular, ghost factors have p values close to zero. 

We depict the connection of j^f\i) with sm{t) via p(N,£) in Figure I5T31 for the number 
to be factorized N = 559 = 13-43 and the truncation parameter M = 2. The upper-left plot 
represents the master curve (t) | (blue line) given by the absolute value of the normalized 
curlicue sum (15. 7p . The function |sm(t)| is even with respect to r, since 

sm(-t) = s* m (t). 

Hence, it depends only on the absolute value of r. Moreover, we note two characteristic 
domains of |sm(t)|: W the function starts at unity for r = and decays for increasing 
r. This central peak around r = is the origin of the ghost factors, (ii) After this initial 
decay oscillations set in whose amplitudes seem to be bound. Indeed, in the Appendix [E] we 
show that in the limit of large M the absolute value of the normalized curlicue sum |sj\/(r)| 
evaluated at non-zero rational r is bounded from above by 1/ \/2. 

The lower-left plot shows the distribution of the fractional parts p(N,£) given by (15.51) . 
The dots in the upper-left plot arise from the projection of the fractional parts (15. 5p of 
the lower-left plot onto the master curve. Those data points represent the factorization 
interference pattern for iV = 559, as depicted on the right. 

5.3 Upper bound on the truncation by complete sup- 
pression of ghost factors 

Ghost factors emerge from the central peak of the absolute value of the normalized curlicue 
function. Our goal is to suppress these ghost factors by increasing the truncation parameter 
M. For this purpose we display in Figure 15.41 the normalized curlicue sum (15 .7p in the 
neighborhood of r = in its dependence on r and M. Indeed, we find a narrowing of the 
central peak with increasing M. In this way we can suppress the ghost factors below a natural 
threshold. 
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N = 559 = 13 ■ 43 




Figure 5.3: Connection between the truncated Gauss sum (£) and the normalized 
curlicue sum %(r) established by the fractional part p(N,£) of the fraction 2N/£. Here 
the number N to be factorized is N = 559 = 13-43 with the truncation parameter M — 2. 
In the lower-left plot we assign to every value of £ the fractional part p(N, £) defined by (15. 5p 
for the number N = 559 as exemplified by £ = 20 and the green dot. In the upper- left plot 
we show the master curve | S2 (t) | indicated by the red curve and given by the absolute value 
of the normalized curlicue sum (15.71) . This curve is an even function with respect to r and 
attains the values above l/y/2 only in the narrow peak located at r = 0. The factorization 
interference pattern for N = 559 shown in the upper-right corner follows from the dots in 
the upper-left plot in a two step process going through the master curve: from £ we find 
the fractional part p(N, t) which determines through the master curve the signal value as 
indicated by the arrows. 
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sm(t) 



0.5 



Figure 5.4: Absolute value |sm(t")| of the normalized curlicue function in the neighborhood 
of r = in its dependence on the fractional part r and the truncation parameter M. The 
function starts at unity for r = and decays for increasing r. This decay becomes faster as 
we increase M. This behaviour is at the heart of the suppression of the ghost factors. 

As shown in Appendix |E] for non-zero positive rational r = p/q the absolute value of 
the normalized curlicue sum is asymptotically bounded from above by 1/ a/2. Due to the 
connection (15.61) between the normalized curlicue sum sm{t) and the Gauss sum J\^f\l) it 
is natural to use this bound as a natural threshold between factors and non-factors. This 
observation allows us to define the ghost factor properly: ghost factors I of a number N arise 
when the fractional part p(N,£) of 2N/£ leads to a value of the normalized curlicue function 
|sm(p)| in the domain between l/y/2 and unity. 

We determine the truncation parameter M such that we can push the absolute value 
of the Gauss sum for all ghost factors below the natural threshold of l/y/2. Ghost factors 
appear for small values of r. This fact allows us to replace the Gauss sum by an integral 
which leads us to an estimate for the truncation parameter M . 

Indeed, with the substitution u = \f2rm we can approximate the normalized curlicue 



function 



M 




(5.8) 



o 



with the Fresnel integral [ST] 
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Figure 5.5: Comparison between the exact discrete normalized curlicue sum (15. 7p shown by 
diamonds and its approximation (15. 8p by the continuous Fresnel integral depicted by the 
black curve. We display the absolute value |sji/(t)| as a function of the number M of terms 
contributing to the curlicue sum (15. 7p for r = 10~ 3 . The horizontal line marks the threshold 
1 / v2 of the signal and the vertical line indicates the upper bound M (I5.10P required to 
suppress a ghost factor corresponding to r = 10~ 3 . 

familiar from the diffraction from a wedge [117] . We note that in the continuous approxima- 
tion the normalized curlicue function depends only on the product M ■ \/2t. 

In Figure 15.51 we compare the absolute value of the discrete curlicue sum s m (t) and the 
continuous Fresnel integral F(Mv / 2r)/(Mv / 2r) at small value r = 10~ 3 . This approximation 
impressively models the results of the discrete curlicue sum. 

We are looking for the truncation parameter Mo such that for a given fractional part r 
the absolute value of the integral ( 15. 8p is equal to We denote a(£) as the solution of the 
transcendental equation 

\F(a)\ 
a 

In particular, for the natural threshold £ = l/y/2 defining the ghost factors we find the 
numerical value of a(£) ~ 1.318. From the fact that F depends only on the product of 
M ■ y/2r it follows that 

a(0 = Mqv 7 ^. (5.9) 
For the factorization of the number N the argument £ is varied within the interval [1, y/~N]. 
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Consequently, the minimal fractional part 

Pmin(A0 



N 



arises from the ratio 2N/£ when the denominator takes on its maximum value £ = yN. 
Finally, we arrive at 

M « - Q ® « ^±VN. (5.10) 
V^~(iV) 2 

Hence, Mo represents an upper bound for the number of terms in the truncated Gauss sum 
(15. 3p required to push all non-factors below the threshold of £. In particular, we find that 
to suppress all ghost factors below the natural threshold £ = l/y/2 we need M ~ 0.659 v^iV 
terms in the truncated Gauss sum. However, we point out that the power-law (I5.10p arises 
from the fact that we use quadratic phases and will be unchanged by relaxing the threshold 
value £, as the change of this threshold will only change the pref actor a(£). 

We conclude this section by noting that the scaling law rests on approximating the normal- 
ized curlicue sum by the Fresnel integral. In Appendix [F] we analyze the range of applicability 
of the Fresnel integral approximation (15. 8p and find that our results lie within the validity of 
the approximation. 

5.4 Ghost factor counting function: inevitable scaling 
law 

In the preceding section we have derived a scaling law between the number M of terms of 
the truncated Gauss sum to factor a given number N. This estimate is a sufficient condition 
for the success of the Gauss sum factorization scheme. In the present section we show that 
it is also a necessary condition. In order to illustrate this feature we first choose logarithmic 
truncation M = In iV and show that at the end of our factorization scheme we will be left 
with too many candidate factors, most of them being a ghost factor. Moreover, we show that 
we cannot achieve a more favorable scaling than the fourth-root dependence, (15.101) . even if 
we tolerate a limited number of ghost factors. 

To answer these questions we introduce the ghost factor counting function 



g(N, M) = # lt= 1, . . . , [VN\ with ±= < \A^\£)\ < 1 J 



(5.11) 



which yields the number of data-points with critical values of the signal in the factorization 
interference pattern for a given N and a chosen truncation M. 
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Figure 5.6: A logarithmic dependence of the truncation parameter M on N is not sufficient 
to suppress ghost factors. The ghost factor counting function g(N, M) calculated for 10000 
random odd numbers iV out of the interval [1,2 • 10 10 ] with M = InN grows faster than 
the logarithm of N. The black line describes the general trend given by (I5.14p . We observe 
strong deviations, as exemplified by N = 13064029441 highlighted by the star and discussed 
in Section [5.4.21 



In order to study the behaviour of the ghost factor counting function g(N, M) on a broad 
range of numbers N we show in Figure \51)\ g(N, M = InN) for 10000 random numbers N out 
of the interval [1, 2- 10 10 ]. Here we choose the truncation parameter to depend logarithmically 
M ~ In iV on N. This result shows that the number of ghost factors g(N, M) for M ~ In N 
grows faster than the logarithm of N. Hence the logarithmic truncation M ~ IniV is not 
sufficient for the success of our Gauss sum factorization scheme. We provide an explanation 
for the general trend observable in Figure 15.61 in Section 15.4.11 and discuss the deviations in 
Section 15X21 

In evaluating the number of ghost factors we proceed in two steps. First, we make use of 
the connection ( 15. 6p between the truncated Gauss sum A^ 1 (£) and the normalized curlicue 
sum %(r). As already pointed out in Section I5T21 the ghost factors appear only for r values 
lying in the small interval [—To, To] around zero. The Fresnel integral approximation from 
Section 15.31 allows us to determine the fractional part tq where the normalized curlicue sum 
assumes the value l/y/2. In the second step we relate the number of ghost factors g(N, M) 
to t by a density argument. 

We determine To with the help of the continuous approximation of the curlicue sum. From 
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(j5.9p we obtain 

To = To(M) w 2^ (5 - 12) 
and we thus arrive at the total width 2r ~ a 2 /M 2 of the interval of fractional parts resulting 
in signal values larger than 1 / a/2. 

We relate the number of ghost factors g(N,M) to the width of the interval 2 To via the 
distribution of fractional parts r for a given N. First, we consider a uniform distribution. 
Here we derive an analytical estimation for g(N,M) which explains the general trend in 
Figure 15.61 Second, we discuss the case of numbers N where the distribution of fractional 
parts cannot be approximated as uniform. Finally, we analyze a trade-off between a smaller 
truncation parameter at the expense of more ghost factors. We show that this approach will 
not change the power-law (I5.10p . 

5.4.1 Uniform distribution of fractional parts 

Let us first assume for simplicity that the distribution of the fractional parts r is uniform for 
a given number N. Here the number of ghost factors g(N, M) is directly proportional 

g(N,M) _ 2r 

to the width 2tq of the interval of the fractional parts which lead to ghost factors. 

Recalling the dependence of tq on M (15.121) we conclude that the number of ghost factors 

g ( N) M)^ l -{^) 2 ^N. (5.13) 

depends via an inverse power-law on the truncation parameter M. 

In Figure I5TB1 we already found indications that g(N,M = lniV) grows faster than the 
logarithm of N. Indeed, from (15. 13)) we obtain 




(5.14) 



which implies that g(N,\nN) behaves like y/~N. In Figure [5761 we display a fit according to 
(15.141) . We find that this fit well describes the general trend over a large range of numbers 
N. However, we also observe strong variations around this general trend. The deviations 
indicate that the distribution of fractional parts is not uniform for certain numbers N. We 
analyze such numbers in Section f5. 4. 21 
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5.4.2 Non-uniform distribution of the fractional parts 

In Figure 1531 we find that for certain numbers the actual number of ghost factors g(N,M) 
considerably deviates from our estimation ( 15.14)) . In the following we show that for such 
numbers the distribution of the fractional parts cannot be treated as uniform. 

This unfavorable case occurs when N has only few divisors, but another number N' = 
N + k close to N has a lot of divisors (with \k\ <C N). For example for the number 

N = 13064029441 = 21647 • 603503 



highlighted in Figure 15.61 by the circle we find that 

N' = N - 1 = 2 s • 3 • 5 • 11 • 17 • 23 • 113 

obviously has a lot of divisors. 

Let us consider £' which is a divisor of N' = N + k but not of N. It follows that if £' > 2k 
the fractional part of 2N/£' is equal to 

2k 

p(N,f) = - T . (5.15) 

If we consider a plot of the fractional part p(N, £) of 2N/£ as a function of £ we will find that 
for divisors £' of N' the resulting fractional parts are aligned on the hyperbola (I5.15P and are 
attracted to zero. Hence for N the distribution of fractional parts p(N, £) is not uniform. 

In the factorization interference pattern of N data-points associated with arguments £' 
corresponding to divisors of N' are also aligned on the curve 



SM 



2k 



(5.16) 



As for large values of £' the associated fractional part —2k/£' tends to zero the resulting 
signal values {£') | approaches unity. Hence the divisors of N' become ghost factors of 
N. 

We illustrate this fact in Figure 15.71 where we plot the distribution of the fractional parts 
and the factorization interference pattern for two numbers: N' rich in factors and N = N' — 1 
rich in ghost factors. To emphasize the region of fractional parts which lead to ghost factors 
we use the logarithmic scale. Here we have chosen N' = 13335840 = 2 5 • 3 5 • 5 • 7 3 which 
obviously has a lot of divisors, as depicted on the upper-left plot by the straight line of 
black diamonds. In the factorization interference pattern shown on the right these divisors 
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correspond to a straight line of signals equal to unity. However, the divisors of N' are non- 
factors of N = iV' — 1 = 13335839 = 11 -479 -2531. Moreover, they are aligned on a hyperbola 
(I5.15P and attracted to zero as shown in the lower-left plot where we can clearly identify the 
hyperbola of stars. Consequently, in the factorization interference pattern plotted on the 
right this hyperbola of arguments with small fractional parts ( I5.15P translates into the curve 
of ghost factors, as depicted on the right. 

5.4.3 Optimality of the fourth-root law 

In Section I5T31 we have derived the fourth-root law (I5.10p as an upper bound on the truncation 
parameter. We will show that it is also necessary for the success of our factorization scheme. 

The analysis of g(N, M) revealed that it behaves similarly to the inverse power in M 
f)5.13p . The closer the distribution of the fractional parts for a given N the better the 
estimation ( 15 . 1 3[) fits the actual data. 

In Figure 15.81 we present the log- log plot of g(N, M) as a function of the truncation 
parameter M for three characteristic examples. First, for the number N = 13335769 which 
has the fractional parts p(N,£) of 2N/£ distributed almost uniformly we find that scaling 
~ M~ 2 predicted by ( 15 . 1 3[) is obeyed. For N = 13335839 we find strong deviations for larger 
values of M due to the fact that the actual distribution of fractional parts is not uniform. 
Finally, for N = 13335840 the ghost factor counting function g(N,M) decays even faster 
than the estimation ( 15 . 1 3[) predicts. Nevertheless, in all three cases the number of ghost 
factors drops down rapidly in the beginning. 

The inverse power-law ( 15 . 1 3|) suggests an alternative truncation of the Gauss sum when 
we tolerate a limited number of ghost factors, say K. Indeed, the power-law reduces the 
number of ghost factors considerably for small values of M. On the other hand, it has a long 
tail, which implies that we have to include many more terms in the Gauss sum in order to 
discriminate the last few ghost factors. However, this approach will not change the power 
law dependence of M, (I5.10p . as the equation ( 15 . 1 3|) yields that 



terms are required to achieve this goal. Let us point out that this results holds if we can 
approximate the distribution of the fractional parts by uniform distribution. However, as we 
have seen in Figure 15.81 if this simplification is not feasible such Mk might be even greater. 
Therefore we cannot achieve a better scaling on N than \fN \ even if we tolerate a limited 
number of ghost factors. 
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N' = 13335840 = 2 5 ■ 3 5 • 5 • 7 3 




I 

N = 13335839 = 11 ■ 479 • 2531 




t £ 

Figure 5.7: Emergence of ghost factors of N from factors of N'. We display the distributions 
of the fractional parts (left column) and the factorization interference patterns (right column) 
for the numbers N' = 13335840 = 2 5 ■ 3 5 ■ 5 • 7 3 which is rich in factors and N = N' — 1 = 
13335839 = 11 ■ 479 • 2531 which is rich in ghost factors. To emphasize the region of fractional 
parts which lead to ghost factors we use a logarithmic scale for \p\ on the vertical axes. The 
number N' has a lot of divisors, as depicted on the upper-left plot by the straight line of 
black diamonds. In the factorization interference pattern shown on the right these divisors 
correspond to a straight line of signals equal to unity. However, the divisors of N' are non- 
factors for N = N' — 1. Moreover, they are aligned on a hyperbola (I5.15P and attracted 
to zero as shown in the lower-left plot where we can clearly identify the hyperbola of stars. 
Consequently, in the factorization interference pattern shown on the right this hyperbola 
translates into the curve of ghost factors. 
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Figure 5.8: The number g(N, M) of ghost factors expressed by the ghost factor counting 
function (15. lip as a function of the truncation parameter M for three characteristic examples. 
We use a log-log plot to bring out the scaling of g(N,M) with M. For the number N = 
13335769 the scaling g(N, M) ~ M~ 2 predicted by (I5.13P with the help of the Fresnel integral 
is satisfied. In contrast for N = 13335839 which is rich in ghost factors we see a strong 
deviation. Finally, for iV = 13335840 which is poor in ghost factors due to the fact it 
has many divisors the ghost factor counting function g(N,M) decays even faster than the 
estimation ( 15.1 3p predicts. 
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We conclude that the scaling M ~ yN of the upper limit of the Gauss sum A N provides 
both sufficient and necessary condition for the success of our factorization scheme. Using 
M terms in the Gauss sum we can suppress all ghost factors for any number N. From the 
relation ( 15. 4 p we see that we need to add 

terms for the success of the factorization scheme based on the truncated Gauss sum. In 
comparison with the value of 1Z ~ N of terms required for the complete Gauss sum (15. 2p 
we have gained a factor of fourth-root. We emphasize that we cannot reduce the amount of 
resources further. 

5.5 Conclusions 

We have analyzed the conditions required for the success of the factorization algorithm based 
on the truncated Gauss sums. Four distinct classes of candidate factors i with respect to 
the number to be factorized N have been identified. In particular, with the help of the 
normalized curlicue sum we have found a simple criterion for the most problematic class 
of ghost factors. The natural threshold of the signal value of the Gauss sum which can be 
employed to discriminate factors from non-factors was identified. We have derived the scaling 
law M ~ yfN for the upper limit of the Gauss sum which guarantees that all ghost factors 
are suppressed, i.e. the signal values for all non-factors lie below the natural threshold. 
Unfortunately, we cannot achieve a more favorable scaling even if we change the threshold 
value or tolerate a limited amount of non-factors. 

However, a generalization of Gauss sums to sums with phases of the form m? with 2 < j 
might offer a way out of the fourth-root scaling law. Indeed, such a naive approach suggests 
the scaling law M ~ 2 i/~N. For an exponential phase dependence m m we would finally 
achieve a logarithmic scaling law. However, these new phases bring in new thresholds and a 
more detailed analysis is needed. The answer to these questions is presented in the following 
Chapter |6j 

Moreover, the analysis of the non-uniform distribution of the fractional parts provides us 
with a new perspective on the ghost factors. So far we have treated them as problematic 
trial factors which might spoil the identification of factors from the factorization interference 
pattern. However, the fact that the ghost factors of N are factors of numbers close to iV 
offers an interesting possibility - by factorizing N we can find candidate factors of numbers 
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N = 32183113 = 613 • 52501 




Figure 5.9: Factors of iV± k obtained from the ghost factors of the factorization interference 
pattern of N = 32183113 = 613 • 52501 with the truncation parameter M = 17 In N. Such 
a choice of M is clearly not sufficient to suppress all ghost factors. However, the remaining 
ghost factors can be fitted to the curves 7^ (£) for k = 1, . . . , 5. Hence, we can identify 
candidate factors of numbers close to N, in our case up to N ± 5. 

close to N. Indeed, as we have found in (15.16P the factors of N ± k align on the curve 7^. M ' ) (I) 
in the factorization interference pattern of N. Hence, if we identify the data points lying on 
these curves we find candidate factors of N ± k. However, to take advantage of this positive 
aspect of ghost factors we need a very good resolution of the experimental signal data. 

We illustrate this feature in Figure 15.91 on the factorization interference pattern of N = 
32183113 = 613 • 52501. Here we have chosen the truncation parameter according to M ~ 
In N 17 which leads to an interference pattern with several ghost factors. However, we can 
clearly fit the ghost factors to curves 7^ (£) for k = 1, . . . , 5. Hence, by factorizing iV we 
also find candidate factors of iV ± k with k = 1, . . . , 5. 
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Chapter 6 

Factorization with Exponential sums 



Introduction 

In the present Chapter we extend the idea of factorization with the help of Gauss sums 
by considering exponential sums. Here the phase is proportional to m? where m is the 
summation index and j is an integer. We show that in such a case the truncation depends 
on the inverse of this function, i.e. M ~ %/N. Hence, we can save experimental resources 
by employing rapidly increasing phase functions. The extreme limit of an exponential sum 
where the phase varies exponentially with the summation index, i.e. m m , should then be the 
optimal choice. We briefly address this case and demonstrate by a numerical analysis that 
the truncation parameter depends only logarithmically on the number to be factored. 

It is interesting to note that recently an experiment |104] based on NMR has used an 
exponential sum with j = 5 to factor a 17-digit number consisting of two prime factors of the 
same order. In this experiment 7r-pulses |118] drive a two-level atom. By choosing the phases 
of the pulses appropriately we can achieve a situation in which the resulting polarization is 
determined by a truncated exponential sum with a particular choice of j. Moreover, even 
the extreme case of an exponential phase m m can be realized in this way. 

We introduce exponential sums in Section [6TT1 and show that they allow us to discriminate 
between factors and non-factors. In particular, we demonstrate by a numerical example that 
phases which increase as m 3 suppress ghost factors more effectively than Gauss sums which 
have phases proportional to m 2 . This feature is our motivation to study the factorization 
properties of exponential sums. In Section l6\2l we have shown that for truncated Gauss sums 
the influence of the truncation parameter M depends crucially on the choice of trial factors. 
We have identified four classes: (i) factors, which are not influenced by M, (ii) threshold trial 
factors, which are also independent of M, (iii) typical non-factors, which decay very quickly, 



117 



and (iv) ghost factors, which decay slowly. In Section 16.21 we perform a similar analysis for 
exponential sums. The numerical calculations of Section 16.11 are confirmed in Section 16.31 
by an analytic argument. We show that the number of terms which have to be summed in 
order to suppress the signal of all ghost factors depends on the 2j-th root of the number 
to be factored. For all exponential sums except the Fourier sum there exist non-factors for 
which the signal cannot be suppressed below certain thresholds by further increasing the 
truncation parameter. The values of these thresholds are determined by the power j and can 
be close to the maximal signal of unity corresponding to a factor. In such a case we cannot 
achieve a sufficient contrast between the signals of factors and non-factors. We discuss the 
restrictions imposed by this fact on our factorization scheme in Section 16.41 Our analysis 
indicates that rapidly increasing phases suppress ghost factors most effectively. This feature 
suggests to consider the extreme case with the phase m m . We briefly address this case in 
Section 16.51 where we present numerical simulations indicating that the resources scale only 
logarithmically. However, in contrast to sums involving a fixed exponent, we no longer have 
the tools of number theory at hand to prove perfect discrimination of factors from non-factors. 
Nevertheless, in the Appendix Owe demonstrate that the sum actually discriminates factors 
from non-factors. We summarize our results in the conclusions of Section 16.61 



6.1 Factorization with exponential sums 

For our purpose to factorize numbers we use truncated and normalized exponential sums of 
the type 



A 



(M,j) 
N 



M + 



1 - 



exp 



m=0 



2ni m? 



(6.i; 



where the phases are determined by the integer power j. Here N is the number to be 
factored and £ is a trial factor which scans through all integers between 1 and J • I n 
the experiments performed so far the upper bound M in the sum is equal to the number of 
pulses applied. 

In the case of j = 1 the exponential sum reduces to a Fourier sum. For j = 2 we find the 
truncated Gauss sum 



(M,2), 



M + 



1 M 

— Y 
+ i ^ 



exp 



m=0 



2 N' 
2irim — 



(6.2) 
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In the case of j = 3 the sum 




M+l 



1 



M |- 

y exp 2ni m' 




(6.3) 



m=0 



is the truncated version of the Kummer sum named after the mathematician Ernst Kummer 
(1810-1893). 

The capability of the exponential sums, ( 16. ip . to factor numbers stems from the fact that 
for an integer factor q of iV with iV = q ■ r all phases in A^f'^ are integer multiples of 2tt. 
Consequently, the terms add up constructively and yield A^f (q) = 1. When £ is not a 
factor the phases oscillate with m and the signal |^4^ f '^(£)| takes on small values. In order 
to factor a number N we analyze for arguments i out of the interval [1, sTN\. We 

refer to the graphical representation of the signal data as factorization interference pattern. 

In Figure lBTTl we show the factorization interference patterns of the number iV = 6172015 = 
5 • 379 • 3257 resulting from the Gauss sum (left) and from the Kummer sum (right) for the 
choice of the truncation parameter M = 15 ~ IniV. In both cases the factors of N lead 
to the maximal signal of unity depicted by black diamonds. In contrast for most of the 
non-factors the signal represented by gray dots is well suppressed. However, for the Gauss 
sum there appear some non-factors, the so-called ghost factors, where the signal indicated by 
black stars is still close to that of a factor. We recognize that the corresponding factorization 
pattern resulting from the Kummer sum does not display any ghost factors. The origin of 
this positive feature lies in the fact that the cubic phase of the Kummer sum shows a stronger 
increase than the quadratic variation of the Gauss sum. 

6.2 Classification of trial factors 

In the preceding section we have shown using numerical examples that the influence of the 
truncation parameter of the exponential sums depends crucially on the choice of the trial 
factors. In the present section we analyze this feature in more detail and identify four classes 
of trial factors. 

For this purpose we start from the decomposition of the fraction N/£ into an integer k 
and the fractional part 



with \p\ < 1/2. Indeed, the integer part contributes only as the multiplication by unity in 
(16. ip and we find 



p(N,£) = --k 
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Figure 6.1: Factorization interference patterns of the number N = 6172015 = 5 ■ 379 ■ 3257 
resulting from the Gauss sum (upper plot) and the Kummer sum (lower plot). Here we 
have chosen the truncation parameter M ~ IniV ps 15. The factors of N, depicted by black 
diamonds, correspond to the signal value of unity. For most of the non-factors, depicted 
by gray dots, the signal value is well suppressed. However, in the case of Gauss sum we 
note that for a few non-factors, depicted by stars, the signal is close to that of a factor. 
Since such arguments can be misinterpreted as factors of iV we call them ghost factors. The 
presence of ghost factors in the factorization interference pattern indicates that the choice of 
the truncation parameter M ~ In N is not sufficient for the Gauss sum. However, the cubic 
phases in the Kummer sum grow faster than the quadratic phases in the Gauss sum. As a 
result, the truncation parameter M = 15 is sufficient to suppress all ghost factor. Moreover, 
some trial factors result in a threshold value of the signal depicted by black triangles which 
cannot be suppressed by further increasing the truncation parameter M. In the case of the 
Gauss sum the threshold is l/y/2 whereas for the Kummer sum it has the value 0.844. 
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where we have introduced the sum 

1 M 

m=0 

This elementary analysis allows us to identify four classes of the fractional part. Indeed, 
we find in complete analogy to the Gauss sums [VI] : (i) for p(N,£) = the trial factor £ 
is a factor of N, (ii) for \p(N,£)\ = tj the trial factor £ results in a threshold value Tj of 
the exponential sum, where the values of tj and Tj are determined by the power j, (in) for 
p(N, £) appropriately away from the origin the trial factor £ is a typical non-factor of N, (iv) 
for p(N, £) ~ the trial factor £ is a ghost factor of N. 

We illustrate the different dependence of representatives of these classes on the truncation 
parameter M in Figure 16721 using the example of the truncated Kummer sum (16. 3p . We find 
signals which are independent of M and equal to unity. They indicate factors. Moreover, we 
note, a rapid suppression of the signal for a typical non-factor. However, for a ghost factor 
the signal is close to that of a factor and we have to include more terms in the sum (16. 3p in 
order to suppress it. Moreover, we find that for certain trial factors £ the signal levels off at 
a non-zero threshold value and thus cannot be reduced at all. 

6.3 Scaling law of the truncation parameter 

In Section 16.11 we have shown that the ghost factors spoil the discrimination of factors from 
non-factors. Fortunately, we can suppress the signal of a ghost factor by increasing the 
truncation parameter M. In this context the truncated Gauss sums were analyzed in [Vlj 
and it was shown that one needs M ~ y/~N terms in the sum in order to suppress the signal 
of all ghost factors considerably. We derive the corresponding scaling law Mj ~ 2 i/~N of an 
exponential sum A^' J) . In [YD the upper bound for the truncated Gauss sum (16.21) was 
obtained by approximating the Gauss sum by the Fresnel integral. We perform a similar 
analysis for the exponential sums. 

Since ghost factors result from small values of the fractional part p = N/ £ — k we replace 
the exponential sum by an integral, i.e. 

M 

A^\£) = S] M \p) » 1 j e^'dm. 

o 

This approximation is justified by the van der Corput method [119] approximating sums by 
sums of shifted integrals. 
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N = 6172015 = 5 • 379 • 3257 
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Figure 6.2: Four classes of trial factors I illustrated by the dependence of the Kummer sum 
CO I on the truncation parameter M. In order to compare with Figure I6~T1 where M = 
15 as indicated by a vertical dashed line we have chosen again N = 6172015 = 5 • 379 • 3257. 
For factors of N, such as £ = 5 depicted by black diamonds, the signal is constant and 
equals to unity. For typical non-factors, such as £ = 10 depicted by gray dots, the signal is 
suppressed considerably already for small values of the truncation parameter M. However, 
for ghost factors, such as £ = 2337 depicted by black stars, more terms in the sum (16. 3p are 
needed to suppress the signal. Finally, for certain arguments, such as £ = 45 depicted by 
black triangles, the signal levels at non-vanishing threshold and it is impossible to suppress 
it further by increasing the truncation parameter M. 

With the help of the substitution m? p = w- 7 and dm = du/^/p we find 
where 

Fj{x) = - I r^"'V„ . 



This analysis brings out most clearly that for small fractional parts p the truncation param- 
eter M and p appear in the exponential sum only as the product M ■ -J/p. 

In order to suppress the absolute value 1*4.^ CO | below a given value £ we have to choose 
the upper bound M according to 

M ■ 4fp = a 

where a is the solution of the integral equation 
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which leads us to the relation 

M = a(Op"*- 

This result shows that the smaller the fractional part p(N, £) of the ghost factor I the 
more terms are required. Since the largest trial factor is of the order of y/~N the smallest 
attainable fractional part 



Pmin(N) 




gives an upper bound 

Mj « a(0pj n « q(0 ^iV (6.4) 

on the truncation parameter M. 

Hence, in order to suppress all ghost factors of iV we require an order of V^V terms in the 
exponential sum . We point out that the scaling law (16. 4p is inherent in the exponential 

sum since the change of £ only modifies the pre-factor a(£). 

In Figure [Owe illustrate the behaviour of for iV = 10 6 + 1 and t = 10 3 

resulting in the fractional part p(N,£) = 10~ 3 ~ 1/yN as a function of the truncation 
parameter M. We visualize the effect of the power j on the suppression of \A^^\t)\ by 
presenting three different curves: (?) black dots correspond to the Fourier sum with linear 
phases, (ii) diamonds represent the Gauss sum, and finally (Hi) stars result from the Kummer 
sum with cubic phases. We find that for the Fourier sum the suppression of the signal is 
extremely slow. Indeed, according to the estimate (16.41) we need Mi ~ y/~N w 10 3 terms in 
order to suppress the signal considerably. On the other hand, for the Gauss sum already 
M 2 ~ \fN Rj 32 terms suffice to reduce the signal, in agreement with (16 .4p . Finally, for the 
Kummer sum the decay of the signal is even faster. We find that M 3 ~ \/~N w 10 terms are 
sufficient to suppress the signal, in agreement with ( 16 .4p . 

In order to verify the scaling law (I6.4p for a broad range of N we have calculated numeri- 
cally the truncation parameter Mj needed to suppress all ghost factors of iV below the value 
£. We have chosen N randomly from the interval [10 4 , 10 20 ] and considered ^ = 0.7. In Fig- 
ure EH we present the results for the Fourier sum (black dots), Gauss sum (open diamonds) 
and Kummer sum (stars). To unravel the scaling law we use a logarithmic scale for both 
iV— and M— axes. The numerical results are in excellent agreement with the estimates ( 16. 4 p 
indicated by the dashed lines. 
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Figure 6.3: Decay of the signal l^jv^ (^)| for increasing truncation parameter M exemplified 
by the Fourier (j = 1), Gauss (j = 2) and Kummer (j = 3) sum. Here we have chosen 
iV = 10 6 + 1 and t = 10 3 resulting in the fractional part p(N,£) = 10~ 3 « For 
the Fourier sum (black dots) we find an extremely slow decay of the signal. On the other 
hand, for the Gauss sum (diamonds) already M 2 ~ yN ~ 32 terms are sufficient to suppress 
the signal considerably. This requirement is further reduced for the Kummer sum (stars) 
to M3 ~ \f~N w 10. We find that our numerical results are in good agreement with the 
analytical estimate (j6.4p . 
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Figure 6.4: Number Mj of terms needed to suppress the signal of all ghost factors of N below 
the value 0.7 for the Fourier sum (black dots), Gauss sum (open diamonds) and Kummer 
sum (stars). To unravel the scaling of Mj with N we use a log- log scale. The dashed lines 



follow from the estimate Mj ~ 



v 



N given by (jBTIj) . 
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6.4 Threshold 



An experiment must also take into account the limited measurement accuracy. Thus for the 
success of our factorization scheme we need a good contrast between the signals of factor 
and non-factors, i.e. we require that the signals of all non-factors are suppressed below the 
estimated measurement error. However, due to the existence of the thresholds discussed in 
Section 16721 this suppression might be impossible for certain powers j. In such a case we might 
misinterpret the signal arising from a non-factor as that of a factor. Hence, such exponential 
sums Atf 1 '^ are not suitable for integer factorization. 

Relation (16. 4p shows that the faster the phase grows the less terms in the exponential 
sum are needed in order to suppress the signal of a ghost factor argument I. However, the 
suppression of the signal might be impossible for all arguments £, as we have seen already in 
Figure 16.21 This feature is closely related to the power j determining the phase. 

The absolute value \A^f'^{l)\ depends on how many different roots of unity we find in 
the sum. These roots of unity are given by 



For the Fourier sum we find all q different roots exp (2nim/q) with m — 0, . . . , q — 1 of 
unity. Moreover, since they are distributed symmetrically on the unit circle they cancel each 
other out. Hence, for the Fourier sum we can suppress the signal \A^ | of any non-factor 
I below any given value by extending the summation range M. 

However, for exponential sums Af> j) with powers 2 < j we are not guaranteed to find all 
different roots of unity. Moreover, since j ^ 1 the corresponding roots of unity exp (2 r nim l pj q) 
are not necessarily distributed symmetrically on a unit circle. Hence, they do not cancel 
themselves completely. In such a case the signal \A^^\tj\ has a non-zero limit as M tends 
to infinity. This limit value determines the threshold and depends on how many different 
roots of unity we find in the sum and their distribution on the unit circle. If we find only few 
different roots of unity which are moreover close to each other on the unit circle the signal 
{A^'^ityl attains values close to unity and cannot be suppressed further by increasing the 
truncation parameter M , even though t does not correspond to a factor of N. 




(6.5) 



where p/q is the coprime rational representation of p(N,£). This is equivalent to 

■ N 

m 1 —q = 0, 1, . . . , q — 1 mod q, 
i.e. the terms in the exponential sum \A { ^' j) (£)\ att ain at most q different values. 
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Figure 6.5: The roots of unity contained in the exponential sums A N ' (£) exemplified by 
the Fourier sum (j=l, black dots), the Gauss sum (j=2, open diamonds) and a higher order 
exponential sum (j=6, black stars). Here we have chosen N = 99 and £ = 7 which leads to 
p(N, £) = p/q = 1/7. For the Fourier sum we find all seven different roots of unity. However, 
in the Gauss sum only four different roots of unity appear. This number is further reduced 
to just two different roots of unity in the higher order exponential sum with power j = 6. 

The fewest possible terms in the sum 4f ^ for a non-factor £ occur if j + 1 is the prime 
number q from the rational representation of p(N, £). In such a case we find from the Euler's 
Theorem (see e.g. Chapter 3 in [120J) 

,• f 1 if q is not a divisor of m 
Tn — < 

[0 if q is a divisor of m 



so m? ■ p is either congruent to p or mod q. With the help of the periodicity m? ■ p 
: in + qY ■ p mod q and the relation (16.51) we obtain for M + 1 being a multiple of q 




= i(i +(5 -i )e -;). 



Hence we find for the absolute value squared 



A M '%)| 2 = "kl + (*"!) cos(^)) 2 + (q- If sin 2 (^)). 



Substituting q = j + 1 we find for p = 1 the threshold value of the sum 
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\{M,j) 



the threshold will always 



For p > 1 or for more than two different terms in the sum .4.^ 
be smaller. 

To illustrate this we plot in Figure IBTBl the behaviour of the signal (-01 as a function 

of the truncation parameter M. Here we have chosen N = 99 and £ = 7 resulting in 
p(N,£) = p/q = 1/7. Hence, g = 7= l- 6 + l and we find that the signal converges to the 
threshold value Ti(6) w 0.953. 
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Figure 6.6: Emergence of the threshold for the exponential sum ,j) with the power j = 6 
for increasing truncation parameter M. We have chosen N = 99 and i = 7 resulting in 
p(N,£) = p/q = 1/7. The signal converges to the value of Ti(6) ~ 0.953 and cannot be 
suppressed by a further increase of M. 

More generally, for prime denominator q = k ■ j + 1 the sum A^ 1 '^ contains at most k + 1 
different terms. For the case of k = 2 an analogous calculation results in the threshold value 

If ( 2tt 



1 + 2j cos 



2J + 1V \2j + l. 

Obviously, for large powers j the values of T 12 (j) are very close to one. 
The above derived results indicate that the exponential sums 

4T J) with powers j larger 

than two can be used for integer factorization only when the experimental data are sufficiently 
precise. For the Fourier sum the signal for any non-factor can be suppressed below any 
given value. However, according to ( 16. 4 j) we have to include a number of terms of the 
order of the square-root of N to achieve this suppression. The quadratic Gauss sum of 
(16. 2p provides a reasonable compromise between the number of terms needed and the non- 
factor discrimination. The gap between the signal of a factor and the greatest threshold is 
approximately 30% which should be sufficient for the experimental realization. The number 
of terms in the sum needed is according to |VIj reduced to the fourth-root of N. 
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6.5 Factorization with an exponential phase 



One way to improve the scaling law might be offered by an exponential sum where the phase 
is not governed by a polynomial as in (16. ip but by an exponential function. This idea leads 
to the sum 

, N 



m=0 



We present a numerical analysis which confirms a logarithmic scaling law. In Section 



16.31 we have found that the number of Mj terms needed to suppress all ghost factors for the 
|(MJ) scales like Mj ~ 2j/ 



exponential sum A N 



N, i.e. Mj is determined by the inverse function 
of the phase evaluated at \/N. This feature arises from the fact that the rising exponent 
prevents the function from accumulating values near unity for small arguments m, as we 
illustrate in Figure IBT71 



TT 








Figure 6.7: Distribution of the roots e 2mm2 p/i (dots) and e 2mmm p/i (stars) of unity for 
quadratic and exponential phase, respectively. Here we have chosen p = 1 and q = 10 4 . 
Since the fraction p/q is small we observe an accumulation of the roots for small values of m 
in the case of the quadratic phase. 



This result suggests that for the exponential sum E^ 1 already a logarithmic number 
of terms M exp ~ In y/N should be sufficient to eliminate all ghost factors. Moreover, our 
numerical analysis summarized in Figure 16.81 indicates that the largest threshold for Sp M 



occurs around the value 0.5. Hence, we can achieve perfect discrimination of factors from 
non-factors. 
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10 10 1 20 10™ 

N 

Figure 6.8: Number M exp of terms needed to suppress the signal | of all non-factors of 
iV below the value 0.7. To unravel the scaling of M we use a logarithmic scale for N. The 
gray line represents the estimate M ~ In y/N. The plot indicates that already an order of 
In y/N terms in the exponential sum £^ is sufficient to find all factors of N. 

However, in contrast to sums involving a fixed exponent, we no longer have the tools of 
number theory at hand to prove perfect discrimination of factors from non-factors. Moreover, 
since the derivative of m m grows faster then m m itself, standard techniques to approximate 
these exponential sums by integrals cannot be applied. Nevertheless, in the Appendix [G] we 
demonstrate that it is still possible to show that the sum actually discriminates factors from 
non-factors by methods of elementary number theory (see |120] for example). 

6.6 Conclusions 

In the present Chapter we have extended the idea of factorization with Gauss sums to expo- 
nential sums where the phase is governed by a power j of the summation index. These sums 
are also capable of non-factor discrimination in complete analogy to Gauss sums. However, 
the truncation parameter Mj needed to achieve a significant suppression of ghost factors of 
the number N scales like Mj ~ 2 y/N . Hence, we can save experimental resources by em- 
ploying exponential sums with large powers j. On the other hand the gap between the signal 
of a factor and the greatest threshold value shrinks as j grows. Therefore, exponential sums 
with large values of j can be used for integer factorization only if the expected imperfections 
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in the experiment are smaller than this gap. 

We have also presented numerical simulations of factoring numbers using an exponential 
sum with exponentially increasing phases. Here the resources scale only logarithmically. 
Moreover, our results indicate that the gap survives. 

Our results also show a connection to two recent experiments |1041 1108] which factored a 
13-digit and a 17-digit number using a Monte-Carlo sampling technique of a complete Gauss 
sum. This method accepts a small fraction of ghost factors and achieves a logarithmic scaling 
very much in the spirit of the exponential phase. 

It is interesting to compare and contrast these two approaches. Ghost factors arise from 
the addition of neighbouring phase factors which only deviate slightly from each other. How- 
ever, when many terms are added the phase factors are distributed homogeneously on the 
unit circle. The Monte-Carlo technique does not add up consecutive terms but tries to collect 
those terms which almost cancel each other. On the other hand, the exponential phase guar- 
antees that neighbouring phase factors deviate significantly from each other and no ghost 
factors can arise. This feature leads to the logarithmic scaling. 
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Conclusions 



Various schemes for factorization of numbers based on exponential sums have been developed 
recently. Their relative simplicity when compared to the celebrated Shor's algorithm results 
in several advantages for the experimental realizations. First of all, exponential sums can 
be easily implemented in various physical systems. Moreover, thanks to the sufficiently long 
coherence times larger numbers can be factorized. 

In the previous Chapters we presented the necessary conditions for the success of the 
factorization schemes based on exponential sums. We found that the number of terms in 
the sum, which directly translates to the number of pulses in the experiment, needed for 
the suppression of all ghost factors is determined by the inverse of the function determining 
the growth of the phase evaluated at the square-root of the number to be factorized. The 
exponential sums with rapidly growing phases are therefore more suitable for the suppression 
of ghost factors. On the other hand, the non-factors resulting in the threshold signal values 
become a significant problem. In general, with the faster growing phase the thresholds 
appear closer to the maximal signal of unity corresponding to the factor. Hence, for the 
successful physical implementation of the exponential sums algorithm for factorization of 
numbers one has to guarantee sufficient resolution of the measured signal. The quadratic 
Gauss sum analyzed in Chapter provides a reasonable compromise between the number of 
terms needed and the non-factor discrimination. Most of the experiments performed to date 
benefited from this fact. 

Needless to say, the simplicity of the analyzed schemes follows from the fact that they 
do not employ entanglement which is the key for the exponential speed-up of the Shor's 
algorithm over the classical ones. Indeed, factorization of numbers based on exponential 
sums relies only on interference. The resources scale exponentially like in the case of all 
known classical algorithms. To improve this scaling law by involving entanglement is our 
next goal. 
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Appendix E 

Determination of Threshold 



In this Appendix we show that for non-zero positive rational r = p/q the absolute value of 
the normalized curlicue sum is asymptotically bound from above by 1/ a/2. This property 
follows immediately from jlllj . Indeed, as shown in [1 llj the asymptotic behaviour of the 
curlicue sum 

M 

Cm{t) = ^exp (i-Km 2 r) 

m=0 

for rational r = p/q depends on the product p ■ q. We find that for p ■ q being odd the curlicue 
is bounded. In such a case the absolute value of the normalized curlicue sum %(r) decays 
with increasing M like ~ M~ l . On the other hand for p ■ q being even the curlicue 

is unbounded and its growth can be approximated by 

|C M (r)|^(M + l)(r -r 1 -...-r At _ 1 ) 1/2 

where 

Tj = (l/rj_i) mod 1 if 73-1 ^ (E.l) 

belongs to the j-th step in the repeating curlicue pattern [11 Ij with r = r. Consequently, 
the limit of the absolute value of the normalized curlicue sum is non- vanishing and for large 
M we can approximate |sjvf(r)| by the finite product 



\sm(t)\ « (t ■ Tx • . . . • r M _i 



,1/2 



We illustrate this feature in Figure [Ej] where we show two different curlicues Cm(t). In 
the upper plot we choose r = for which the product p ■ q is odd. In such a case the 



8 



function Cm(t) is periodic in M and the curlicue is bounded. On the other hand for r — nnn1 
the curlicue depicted in the lower plot is unbounded and its ultimate growth is linear. 
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Figure E.l: The behaviour of the curlicue sum CmG°/<?) in dependence on the parity of the 
product p ■ q. In the upper plot we choose r = p/q = 9/10001 for which the product p ■ q is 
odd. We find that the curlicue repeats itself and is bounded. As a second example we choose 
t = p/q = 8/10001 where the product p ■ q is even. In such a case the curlicue expands, as 
depicted in the lower plot. 
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Let us determine the asymptotic bound of |sjy(r)|. The recursion (lE.ip terminates [lllj 
at r M = which implies r^-i = 1/b where b is a natural number. Since Tj < 1 we find the 
estimate 

\sm(t)\ <y/f^=l/Vb. 

The case 6=1 cannot be produced by the recursion formula since all Tj are strictly less 
than one. As a consequence the absolute value of the normalized curlicue sum is 
asymptotically bound from above by 1/ Vz. 
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Appendix F 

Applicability of the Presnel 
approximation 

Let us comment the range of applicability of the continuous approximation (15. 8p . The scaling 
law Mo ~ \/~N connecting the number to be factored with the truncation parameter Mq 
necessary to push all ghost factors below the threshold 1 / \/2 relies on the approximation of 
the normalized curlicue function by the Fresnel integral. For large values of N the scaling 
law requires large values of Mq. However, for large M the continuous approximation might 
not hold any more. 

For the continuous approximation to hold the phase difference 

7T ((m + l) 2 - m 2 ) r = 7r(2m + l)r 

of two successive terms in the sum (15. 7p should at most be of the order of ix. Together with 
the fact that the maximal phase difference appears for m = M we arrive at the inequality 

t(2M + 1) < I. 

Indeed, this condition is violated for sufficiently large M. 

When we recall that for a given N the smallest fractional part is r m i n = 1/ y/N we arrive 
at the constraint 




on the maximal value M c of the truncation parameter for a given N. Thus M c ~ \fN provides 
an upper bound on the validity of the Fresnel approximation, (15. 8p . Since Mq ~ \fN the 
Fresnel approximation is valid. 
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Appendix G 

Discrimination Property for Variable 
Exponents 



In this Appendix we prove that the exponential sums with exponential phase allows us to 
distinguish factors from non-factors of a given number. The discrimination property of the 
exponential sums with a fixed exponent rests on the fact that only for integer values of I 
which are factors of N, the sum can take the value unity There is a number theoretical 
argument supporting this fact, as long as the exponent j in the sum (16. ip is fixed. This 
feature comes from the distribution of the values exp(27rim : 'y) on the unit circle. For fixed 
j, it is impossible to hit the same point twice as m increases provided we use a truncation 
parameter M below yN . However, for a variable power m m that is an exponential phase, 
this non-recurrence property is not obvious. In this case we need to prove the discrimination 
property explicitly. 

The value exp(27T2m m y) depends on the fractional part of m m y only. We hit the same 
point twice for different values m and n if and only if 



where k is an integer. 

As in (16.51) we make use of the coprime rational representation of p(N,£) = p/q and find 
that the phase factor 



is a g-th root of unity. In particular, it is the (pm m )-th one if we enumerate them counter- 
clockwise starting from the zeroth root 1 = exp(27ri^). Note that c-th and d-th roots coincide 
if and only if q is a divisor of c — d. 



m m__ n n_ = k 



(G.l) 
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So the discrimination property depends on the fact, that there are values c and d such 
that 

q is not a divisor of pc c — pd d . 

The discrimination threshold does not depend only on the number of such pairs, but also on 
the position of the corresponding roots of unity. Opposite roots of unity eliminate themselves 
in the sum, so the worst case occurs if these roots accumulate on the same position. 

We consider two cases: for large q, the first numbers in the sequence pm m will be below q, 
so any pair chosen from the beginning of the sequence cannot fulfill the recurrence condition 
(IG.ip . so they correspond to pairwise distinct roots. As a consequence, the absolute value of 
the sum cannot assume the value unity. 

For small q, we show that the p-th root exp(2ni-) and its conjugate exp(— 27u|) appear 
in the sum, which leads to the elimination of their imaginary parts. According to Euler's 
Theorem |120] there is an even m such that pm m corresponds to the first root exp(2ni^) and 
j = m/2 gives pji, which corresponds to the (— l)-root exp(— 2ni^) = exp(27ri^-). The sum 
of this conjugate pair is a real number strictly below unity. 
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